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Der Spannungszustand in einem aus zwei elastischen Halbscheiben 
bestehenden Verbundkérper 


Von H. Bufler 


1. Einfiihrung. In einer friiheren Mitteilung! wurden aus den Grenzbedingungen des gemisch- 
ten Randwertproblems die Integralgleichungen fiir die Druckspannung p und die Schubbelastung q 
im Verbundgebiet —a< x <a zweier verschiedener elastischer Halbscheiben gleicher Dicke 
(Abb. 1) aufgestellt und gelést. Die unter Voraussetzung eines ebenen Spannungszustandes * 
gewonnene allgemeine Lésung lautet ** danach in komplexer Schreibweise und in bezug auf Abb. 2 


tt \qhl= gt (a D, + x Dy) (-—a<x<a). (1) 
Hierin bedeuten : 
F(x) = p(x) +iqe) = (i =/—1) (2) 
und 
a a t+ Cy 
Ba eae (3) 
sowie 
5 ie ; . 
3 Dy = 2G 6 +i( T4248), D,=C,+71C,, (4) 
_% (l—)Ej—Q—) By /_ Ed 
Ry a EH, EF, i 
Sareea Tae ee (5) 
a. E, Ee ; Cyst Spi 
er ES (62— a) # + ee 


Der Elastizitatsmodul, die Querdehnzahl und der Warmeausdehnungskoeffizient der Scheibe 1 
bzw. 2 sind dabei mit E,, ,, «, bzw. Ey, ¥2, %2 bezeichnet; t ist die Temperaturdifferenz gegeniiber 


Abb. 2. Randspannungen p, q Abb. 3. Die verschiedenen Beanspruchungsarten 
an der unteren Halbscheibe. j des Verbundkérpers. 


ch. 39 (1959) S. 218. aie ‘ 

1 sta d oes ne Nahe der Berithrungszone beider Scheiben wegen der hier 
Querdehnung (senkrecht zur Scheibenebene) eine gewisse Storung, die auf einen 
bt, je geringer die Scheibendicke ist. Bei sehr dunnen Scheiben, die 
ihr Einflu8 vernachlassigt werden. 
n friiheren Arbeit legen wir hier das Koordinatensystem der Abb. 2 — ‘ie 
veckmaBig ist. 

Ca eae i 
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dem wirmespannungsfreien Zustand (wobei konstante Temperatur in beiden Scheiben voraus- 
gesetzt wird) und zahlt positiv bei Erwarmung; gp bedeutet den Verdrehungswinkel zwischen 
beiden Scheiben, und P baw. Q ist die Kraft pro Scheibendicke, mit der beide Scheiben aufeinander- 
gepreBt baw. gegeneinander geschoben werden; o,, und og sind schlieBlich die konstanten Vor- 
spannungen in x-Richtung. Das Verdrehungsmonemt (pro Scheibendicke) M, die Querkraft Q und 
der Winkel p (bzw. die Konstante C,) sind durch die Beziehung 


7 t2h4) (6) 


2 
COT ep ( 


miteinander verkniipft. 

Gleichung (1) gilt fiir fiinf verschiedene, in meiner friitheren Arbeit diskutierte Belastungsfalle 
(Abb. 3), namlich reinen zentrischen Druck bzw. Zug (Fall I), Schub bei verhinderter Verdrehungs- 
moglichkeit (Fall IT; wegen p = 0 ist dazu das Moment M = — 2 f Q a erforderlich), reiner Schub 
(Fall III; wegen M = 0 stellt sich hier der Winkel gp = : i are “oe : ein), reine Verdre- 
hung (Fall IV), Warme- und Vorspannungen (Fall V). Es wird angenommen, daf sich die beiden 
Halbscheiben auBerhalb des Bereichs — a < x < a nicht beriihren. Die Berandung einer Scheibe 
muB8 dann etwa so verlaufen, wie sie in Abb. | angedeutet ist. 

Trennt man Gleichung (1) in Real- und Imaginarteil, so folgt wegen (2) 


A cos (2 in2—=*) — B sin(@ in 24) 


yat—a? a (V) 


—x : a—x 
_ Gofx8 a. Asin (pin |) 


7 Ja—x8 


ple) = 178g 


a 


B cos (2 In 


a 


q(x) 


A und B sind dabei die spater in (23) festgelegten Abkiirzungen. 

In der vorliegenden Arbeit soll der Spannungszustand im Inneren der beiden Scheiben bestimmt 
werden. Man hat dazu die entsprechende Randwertaufgabe (erstes Randwertproblem fiir die durch 
(1’) gegebene Randbelastung) zu lésen. 


2. Bestimmung der Spannungsfunktion. Zur Ermittlung. der Spannungen in beiden Halb- 
scheiben macht man zweckmabig von der komplexen Methode Gebrauch, die auf G. V. Kolosov und 
N. I. Muskhelishvili! zuriickgeht und deren Anwendung wegen (1) und (2) naheliegt. Danach lassen 
sich die Spannungskomponenten in einer Halbscheibe, deren Randspannungen oder Randverschie- 
bungen vorgegeben sind, durch eine einzige komplexe Spannungsfunktion @(z) der komplexen 
Veradnderlichen z = x + iy ausdriicken. 

Wir betrachten zunachst die untere Halbscheibe entsprechend Abb. 2. Bei gegebener Rand- 
belastung p, q erfiillt die Funktion* 


il Lied =e . 
B,(2) = D(z) = ,* | ee a (7) 
die Randbedingungen 
Cis tT=—q fur = a2 — a? , 
o, =0, GU fir x? >a?. 
Wegen (2) kann man schreiben 
fa 
head! F(u) a 
M0) =a; | Tau, " 


—a 


wobei hier wie im folgenden der Querstrich den jeweils konjugiert komplexen Ausdruck bedeutet. 
Unter Beriicksichtigung von (1) folgt nun 


$0.3 Gt os (7”) 


1 N. I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticit Groni 1953 
* Siehe Gl. (113.2) des oben zitierten Buches, mit unseren Beethro ears 
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mit den Abkiirzungen 


4a ay (eA 
a 
a f ce 
Ji = a du = - aoe — d ut 
Va oats (t 2) he i u\2 fu Zz a 
xv oe a a a 
und 
+a 4 Y tote aie 
a u 
(eee ae Ve u oa 
7 Fase eS MER Ee 


Je—w | \:-(2F 


Wir fiihren folgende Substitutionen ein: 


u eu* — ] ie 
at] 78 3 (8) 
bs ez* —] oe 
ee = oo (°) 
2 
Wegen i) <= List der Hauptwert von u* reell; z* ist im allgemeinen komplex, sein Hauptwert wird 
nur reell fiir z =x und el < 1, also langs der Beriihrungslinie. Aus dem Ausdruck 
a 
eo Op ey) ee 2ay 


~ (a—xP ty tt @—xpesty?? 


der unmittelbar aus (9) mit z =x -+41y folgt, ersieht man, da fiir Punkte innerhalb der be- 
trachteten (unteren) Halbscheibe wegen y < 0 die Beziehung Im (e?") < 0 gilt und infolgedessen 
fiir den Hauptwert von z* auch die Beziehung 

goa lines <a -2 7E-. 


Mit (8) und (9) erhalt man 


“ (p+ iA) 
el ro du * (t < Imz* < 22) (10) 
2a 


und 


(S+i0)u* 


+ oo 
te oe 
ye if eres ae (11) 


— oo 


Zur Bestimmung des Integrales J, dient die von H. G. Séhngen [dort Gl. (19)] abgeleitete Formel 


+ oo 
Teg eRe cet Ses par ee I de 1 (u* reell; 0< Res <n; (12) 
a (eu* +1)” n—1}sinas n ganzzahlig und positiv). 
Sie liefert das Ergebnis 
— 4 ll’ 
Ji Gone; (11’) 


Das Integral J, kann nach der Residuenmethode ausgewertet werden. 


1 H. G. Séhngen, Math. Z. 60 (1954) S. 31. 
: ; 16* 
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Wir berechnen zunichst das etwas allgemeinere Integral 


+00 ; 
1 | aaeile das 
I eu* — @ 
unter der Voraussetzung, daB 0 < Imz* < 2” und 0< Res< 1. Betrachtet man hierin vor- 
itbergehend u* als komplexe Veranderliche, so besitzt der Integrand f(u*) = g(u*)/h(u*) im Gebiet 
Im u* > 0 Pole erster Ordnung bei 


U=2* +2imn (n::0; 122) 


Nach dem Residuensatz gilt dann fiir das sich tiber die gesamte Halbebene Im u* = 0 erstreckende 
Umlaufintegral 


fee) z(u*) . (oe) z 
yee * * 97; a) tt ge e—2insn — —s(2"—i7) 
p f(u*) du ae > veal =f i? » sins ~ 
n= r PS ycee 


Man kann den Integrationsweg folgendermafen aufteilen: 


f flr) du* =f flut) du + ff(Rey) Rie dp. 


— oo 
reelle Achse R- 
Halbkreisbogen 


Da der Halbkreisbogenanteil fiir R—> oo und 0< Res < 1 verschwindet, folgt sofort 


-+- co 


1 * e(1—s) u* rae ae e—s (z* —iz) Qe I % 5 gu 0 R 1 
| eur — et sin 7 Ss” ( ch Re ded hi < Res< 1). (13) 
Aus (10) mit (13) bekommt man nun 
tity 1 res core 5 
LISS Seat er cau rea . (10’) 


Wir setzen jetzt (10°) und (11’) in (7”) ein. Die fiir die untere Halbscheibe giltige Spannungs- 
funktion lautet dann bei Beachtung von (9) mit @®(z) = D*(z*) 


D*(z*) = oe i et b @* — ia), (P, Coj = + D, Sin = ae D, (14) 


oder, in 2 ausgedriickt, 


ee 
aap = = \2-+ a = ' 
D(z) = eG ig (a D, + 2D,) poe + D, | ; (14’) 
Im weiteren wird der Differentialquotient ®’(z) = - @(z) benétigt. ZweckmaBig geht man 


dazu von (14) aus und bildet zunachst 


d®*(z* LOCC RDO Dian tag ee lee ie. Dy ati oe Ve gt 
a 8 a= ( A ip D,) Gin > + (3 + ip Di) Gof = (15) 


sowie 


d Gr a per eiurae ei B (2* —iz) 2*[/D, pute arias Meee 
AI gel ears or =| + ip D,) cin +(2 + ip D,) oo 3]. (16) 


Bei Beriicksichtigung von (9) geht (16) tiber in 


ey ea 
P (ES : i ae jp 
Oe) ae (22 — a®)3/2 (2 ip D,)2+ (2 +ip Dia . 


(16’) 


3. Ermittlung der Spannungskomponenten in der unteren Halbscheibe. Die von der Rand- | 


belastung p, q herrithrenden Normalspannungen (6,)9, (@y)o, und die Schubspannung T, in der 
unteren Halbscheibe sind mit der Spannungsfunktion ®, = @ verkniipft durch die Gleichungen* 


()a + (0,)2 = 2 [D(z) + BZ)] + oy, (17) 
()a— (Og + 2 it, = 2 [@ — 2) ©(z) — D(z) — O(z)] — a, (18) 


* Siehe Seite 454 des in FuBnote 1 von Seite 234 zitierten Buches. 
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(in vorstehenden Formeln ist die zum Rand parallele konstante Vorspannung 0, bereits beriick- 


sichtigt). Dabei ist entsprechend der Definition folgender Zusammenhang zwischen @(z), ©(z) 
und @(z) zu beachten: Falls 


D(z) = D(x, y) + i Py,(x, y) , 
ist, dann bedeutet | 


D(z) = D,(x, y) —i D,,(x, y) 


und =a a 


P(z) — D(x, —y) ay: D(x, —y) é Abb. 4. Bipolarkoordinatensystem. 
Fiir das weitere ist es giinstig, das in Abb. 4 skizzierte Bipolarkoordinatensystem einzufihren. 

Hierfiir gilt 

z+ ary eri, s—a=r,e”, (19) 


Damit erhalt man nach geeigneter Zusammenfassung aus (14’) und (16’) 


D(z) + BE) == 2 (a A + Cy) sin (8 In oS 
+ (a B—y C,) cos (2 In 2 ol - ea — = ; (20) 
D(z) + D(z) = cs ae - ( A +tyC,) sin ip In ie +18 (@,— o)| 
+ (aB + iy C,) cos [Bin +18 (m—o)]|—&, (21) 
; To 3 
2 —t]| pln— + —(@, + @,) 
ee gape ares | (C—iD). (22) 


In vorstehenden Formeln wurden als Abkiirzungen folgende reelle GréBen verwendet: 


fn Ce 
Ba? Giese, ; | 
as 7 t (23) 
C=(+4ANG tT +(F +28) 5, 
D=(1+4—) ¢C, ae 


Setzt man nun die Ausdriicke (20), (21) und (22) in (17) und (18) ein und trennt in Real- und 
Imaginarteil, so erhalt man fiir die Spannungen (0,)>2, (o,)2 und T. in der unteren Halbscheibe die 
drei Gleichungen 


a? ef (@2— ,) 


> 3 
Fs eee {¢ an p In “8 + Se + | + Deos ip Int + 5 (0 + o)| 


a 5 Oe SE, @, + Ws 
a m Vr, 2 le ee os. Gene mae Ye 
. 2 eB(o,.—a,) : : T. | O1+ Oe 
(n)e + (on =F |@ A ty G2) sin(Bin + AS) 
a . 
+ (aB—y G) 0s (Bm @ + FS) — 2 C+ oun, (25) 
al 


D ef (®,— @1) 


r Tp 3 
(0y)2— (Ox)2 = Fe Oy ie pe P sin pin + = (@, + o)| — Coos ip em es (@, + on)| 


bere A (sin 5 + Eco 258) + 3 C—O (26) 


7 Yr, Yo 
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Hierin bedeuten 
j : 2 2 Ta 
ie z Sin B (w, — @)) IC. sin (3 In =) — C, cos ( In | 


1 


“ 
Lal 


4+ Cof B (@,— @) 4 sin (8 In a + Bos (6 In | 


F =~ Goj 6 (o. — 0) IG sin (Bn i) + Cy 00s ( In 2) 


9 | “3 
+ Gin B (w.—@,) = B sin (8 In =| +A cos (é In 7) : 


4, Die Spannungskomponenten in der oberen Halbscheibe. Aus der Lésung fiir den Spannungs- 
zustand in der unteren Halbscheibe 1aBt sich leicht diejenige fiir die obere Halbscheibe her- 
leiten. In bezug auf das (&,7)-Koordinatensystem (Abb. 1) geben namlich die Formeln (24), 
(25) und (26) die Spannungen fiir die obere Halbscheibe, wobei jetzt z = € + 17 bedeutet. Aller- 
dings ist, da beide Ebenen ihre Rollen vertauschen, wegen (5) zu ersetzen 


C, durch — C, (daher —/ statt +), 
C, durch — C,, C, durch — C,, Oo durch 69, - 


Will man die Spannungen wieder auf das x-y-Koordinatensystem beziehen, so hat man weiterhin 


zu ersetzen: 


z durch —z (d. h. x durch — x, y durch — y, r, durch rg, r, durch r,, ©, durch @,—2. @, durch @, —2). 


Damit erhalt man die endgiiltigen Formeln 


a2 eB (:— @)) nage Pete a 3 ey eee 3 , 1] 
Gh CELE {c sin G In = -+ 5 (om, 4 0) _ Decos t In as (@, + 2) 
eer = en ge Oe Fos @, + ®s\ (28) 
a Vr, re - : : 
2 ef (@,.—@,) aa s D1 + @s 
(6) + Ch =— = [lad + y Gy) sin (Bint + 2%) 
ue TT, 1 = 
: 2 x 2 2 . 
+ (a B—y C,) cos (8m - Ee w =) 4 z Ceines 29) 
2 2 8 (@,—@y) . 3 T2 3 ] 
(0): — (¢,); = ah TS 8 = ry SP | sim i, In = = i > (@, + @g) — Ccos |B In y a > (@, + @9) I 
a3 = == (F sin“! z ©: 4 Ecos as — = Cy — 091 - (30) 
bee J 


A, B, C, D, E, F sind dabei identisch mit den friither schon verwendeten Abkiirzungen. 

Die Formeln (28), (29), (30) unterscheiden sich von den Formeln (25), (26), (27) lediglich durch 
das umgekehrte Vorzeichen (wenn man von den konstanten Vorspannungen Oo, und Gog absieht). 
Dieser Sachverhalt ist in der Tatsache begriindet, daB zwischen der Spannungsfunktion der oberen 


[,(z)] und der unteren Halbscheibe [@,(z)] die Bezichung 
D,(2) + Dz) = 0 


gilt, welche sich aus den Randbedingungen herleiten laBt*. 


5. Diskussion des Spannungszustandes. Wir wenden uns nun den in Ziff. 1 besprochenen Be- 
lastungsarten I bis V zu. Die in den einzelnen Fallen maBgebenden GréBen C,, C, [mach (5)], 
A, B, C, D [nach (23)] sowie E, F [nach (27)] sind in Tabelle 1 zusammengestellt. Sie vereinfachen 
sich wesentlich fiir 6 = 0 (Tabelle 2). 


* Vgl. hierzu Seite 491 des Buches von FuBnote 1 von Seite 234. 
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Tabelle 1. Die fiir die Beanspruchungsarten I bis V mafgebenden Gréfen C,, C,, A, B, C, D, E, F 


C, C, A B C D 
1? P P 
I 0 0 = eter Peco R lc. 
= g : a a (2 ; 28) a 
II 0 0 0 Zs Gere ee: 
a a a aoa 
4 = a NS 
J eae . x 4B Q| 1—467Q (2 any e 2@ 
PAR S al+4f?a 1+ 4B? a a a aa 
2M 4 ONE ae AME | uM ; 
IV —- = = 
1 + 4 ? a? . a1+ 42a? | 1148? @ 2 | 0 
| ee } x ee c | re i 
V 0 (Ce mach (5) —28G =a ‘ a 4G 
E 
IT | Gojp (2 — 0) sin (Bn 72) 
yee 
II | of (w,—«,) cos (6 n=) Na 
ee me 


ED \atP< inp (o,— o4)e05(B In 7) — Goff (on —en) [48 sin (B In) —0—4p")e05[PIn 72) 
1 a Ty ry | == 


~ SinB (o,— 0) c08 (2 z Sy Me ae [2 sin(p in?) Bc lr “I ane 


J 
| 
te Sin P (w, — w,) sin (2 In 2 + Co} 6 (mw, — o,) ie cos (e In #4 —26 sin ( (a In al} Ge 


F 
Hf | Sin B (w2, — ,) cos (@ In = = 
II —GinB (w,— o,) sin (e In mie 2: 
III | — feo (o—0,) sin(p nj |+ SinB (w,—o,) 40 eae ( nz) 4 1 —462)sin (6 8 Int? fe Oj = 


IV {Fo (@2— @,) sin (em) + Sin B (w@,.— ©) lz cos (am = ale sin (@1m rea 


") 
sl 
1+ 
an = Cof B (w, — a) cos (Bn 72) — Sin B (w, — oy) 2p cos (@int) + sin (ém alike 


Tabelle 2. Die fiir die Beanspruchungsarten I bis V mafgebenden Gréfien C,, C,, A, B,.CeD, E, F, 
speziell fiir B = 0 


Clin A B Coheed E F 
I | 0 | 4 = 6 & = oo | 0 0 
| | | | | 
FETED |. 5-50 | See eee os | we = £ | sae Q 0 
IV | = a | 0 Sethe 0 | — a | 0 0 0 
V | erate tae) iron | SG, | 0 aed = G | 26, 
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Wie in Ziff. 4 gezeigt wurde, sind die Formeln fiir die entsprechenden Spannungen in beiden 
Halbscheiben in bezug auf das Koordinatensystem der Abb. 4 identisch bis auf das Vorzeichen, 
also 


Tol Ta» Xo Vo Ops Way B) = — Ty (Tys Ta» X Vo M1» Da» B) » 
(Gy)o (Tr Tar % Yo Dy> Oo, B) = — (Gy)1 (Tas Far Xs V2 Or» Oo» B) » , (31) 
[(o.)2 — Goa] (Tre To» Xs V2 Oy» Mo» P) = — [(ex)1 — Ooi (Tr, Too + V2 M1» Wa» B). | 


Die Veranderlichen r,, rg, 1, >, x, y sind entsprechend Abb. 4 miteinander verkniipft (wobei fiir 


Oy 
T, (0x)1» (0,), (obere Halbscheibe) stets my ¢ = 0 und fiir 7, (o,)o, (Gy)z (untere Halbscheibe) 
4 yx 


stets w, ¢ <0 gilt). Der Parameter f [nach (3)] beinhaltet die elastischen Konstanten der 


af 
beiden Scheiben und liegt im Intervall 
In 3 In 3 
Fe SPS on 
denn wegen (5) ist 


IU a4 
Rs SE Ge 
7 SQ85- 


Falls Scheibe J (bzw. 2) starr ist im Vergleich zu Scheibe 2 (baw. 1) und Scheibe 2 (bzw. 1) die 
J In 3 
Querdehnzahl Null besitzt, gilt Cj = + (bzw. —) + und daher 6 = + (bzw. —) = 


Qn? 
dagegen beide Scheiben dieselben elastischen Konstanten, so gilt C) = 0 


besitzen 


ie : 
S und infolgedessen 6 = 0. 
z y Die iibrigen in den Formeln fiir die Spannungen auftretenden 
GréBen sind die resultierenden Krafte P,Q und das resultierende 
iy / ; Moment M sowie der im Fall V (Warme- und Vorspannnngen) be- 
y nétigte Ausdruck C,. 
p Ks soll jetzt der Zusammenhang zwischen symmetrisch zur x-Achse 
/, 
ips sos kas gelegenen Punkten P und P’ untersucht werden (Abb. 5). Bedeuten 
legene Punkte P und P’. ' 
Rese =, O=—m, yay, f'=—f, (2) 


so erhalt man aus (24) bis (30) und mit Hilfe der Tabelle 1 


Ty = Ty (Tas Too so gs Vas Ps 


T: = — fT (Ty Te, ¥, Y's Wy, Wo, B’) , 
(oy)1 os (oy) (Ty. Tas XY, My, Wo, B) » 
(6y)2 =f (ey) (71, Te. X,Y", W4> Wo, B’) » (33) 


2 
(6), — G1 — = C, = 


2 
(Ox)1 — %1 EEF c,| (Ty, Te. X Yo Wy, Wo, B) 


2 o) eas 
(6,)2 — Go + oe C, = f |(6,)1 — o1 ie Cy] (71, To, XY’, Wy) Mg, B’) « 


Der Faktor f regelt hierbei das Vorzeichen und bedeutet 


a] 


f= i 1 bei den Fallen | tle 


TELE. (34) 


Wir haben somit folgendes Ergebnis, das auch fiir die numerische Auswertung niitzlich ist: Aus den 

Formeln fiir die Spannungen an einem beliebigen Punkt P der oberen Halbscheibe erhalt man 

diejenigen fiir die Spannungen an dem zu P beziiglich der x-Achse gespiegelten Punkt P’, indem 
a ot 4 

man lediglich 6 durch 6’ = — B sowie w, ¢ durch } ws ersetzt; das Vorzeichen der Spannungen 
. . . a. oe 

wird dabei durch (33) und (34) festgelegt. Im Fall = 0 ist daher dem Betrage nach kein Unter- 


schied zwischen den entsprechenden Spannungen in P und P’ mehr vorhanden. 


: 
} 
z 
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Fiir 6 = 0 werden die Formeln besonders einfach. Sie lauten im Fall I 


oh eet: es, 3 
Ls eae (r, r2)3/2 pelea ts (@, + 2) — y cos x (oy + @s)}, 
2 Ge ieee 3 | 
(o,)4 m Wir, [sin 5 (w, + ws) one y sin — 3 (Wy, + Ws) + x cos a (@, + o)| ; (35a) 
Pil a oan ars : , 
(o,),; = — = inn [sn 5 (@, + @) + ae y sin > (w, + m,) + x cos : (a, o)} J 
im Fall If und ITI 
Ostett ESS 
Ui lee eae rag fring (@, + @,) + ae ly sin 5 (@, + w,) + x cos — > (04 ++ | 
2 
Q ¥y was 3 
(o,), =— 7m re |* 8 (w, + @,) — y cos = (a, + ,)), (35b) 
eae act 1 <a ee 3 
6); = = aes cos > (@, + @s) nes x sin 5 (@, + @,) — ¥ cos 5 (wm, 4 | : | 
im Fall [IV 


2 ¥y eS 
Ty ee Gere oO ++ Ws) » 


(oy), aR ft 


(,), == M—— | - 


| 
x sin + (o + @,)— y cos (a, 4 9) ¥ cos (wm, 4 o| » $ (35e) 


ee! L: 3 
« sin — (@, + 2) — y¥ cos 5 (w, 4 | t+ cos 5 (a, 4 oh 
it 


Ae To YT, To 
im Fall V 
Greta if Al art 3 
ce ae ie x sin > (@, + 2) —y cos > = (0 a 0) + 77008 y (1 + on} 5 
C a sibs 
(¢,), = ae (yrs 2” as (@, + @s) ; 
pee (35d) 
a 3 Set) 
G.)5 = —. eS la ly sin 5 sin — x (oy + W,) + x cos > (oy + o,)| — = sin > (@, + on} 
2 
+- Ga C, 7 001 . J 


Die Gleichungen (35a), (35b) und (35c) sind auch als Spezialfall in den Formeln fiir den Spannungs- 
zustand in einer unendlichen Scheibe mit beiderseitiger hyperbolischer AuBenkerbe (namlich bei 
verschwindendem Kerbradius) enthaltent. 

Wie aus den in Ziff. 3 und 4 abgeleiteten Formeln hervorgeht, wachsen die Spannungen bei 
Anndherung an die Eckpunkte (a, 0) und (— a, 0) wegen r,—> 0 bzw. r,—> 0 iiber alle Grenzen. 
Dieser Effekt ist bereits von den Spitzkerben her bekannt!. Auf eine genauere Untersuchung der 
Besonderheiten des Spannungszustandes in unmittelbarer Umgebung der Eckpunkte wird im 
_ nachsten Abschnitt eingegangen. 


6. Auswertung an speziellen Stellen. In Tabelle 3 sind die fiir die Koordinaten maBgebenden 
Relationen zusammengestellt, und zwar unter der Spalte 

a) diejenigen, welche sich auf die x-Achse bezichen, 

b) diejenigen, welche sich auf die y-Achse beziehen, 

c) diejenigen, welche sich auf die Gerade x = a beziehen, 

d) diejenigen, welche sich auf den kleinen, um den Eckpunkt x = a gelegten Kreis vom Halb- 
messer. 9 beziehen (g/a <1). Die Zeichen © und © sollen hierbei andeuten, ob die obere oder 
untere Halbscheibe gemeint ist. 


1 H. Neuber, Kerbspannungslehre, Berlin/Géttingen/Heidelberg 1958. 
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Tabelle 3. Beziehungen fiir die Koordinaten. — a) langs der x-Achse; b) langs der y-Achse; c) langs 


der Geraden x =a; d) entlang des um den Punkt x = a gelegten Kreises vom Radius 0 2 < 1) 
eee 
a | 
oe See tab) ¢) d) 
Wee | > a | 
® 2) | 2) | 2) | @ 2) ® | ® @®; @ 
ss Sn nn ee nn 
x x 0 a | a 
y | 0 | atg @, 2atgo@, | esin@, 
@, | 0 } 0 0 | Wy, . WM, 0 
a A hoes Bos erp +. : : | ears 
| | Pe __, tir ¥z 0 | ie | ar 
Ws Ir | re | 76 —4|—H— Oy] | ak Ws 
- =e . ae 
ae are ae 2.8 | 1 isinnced,f orig eie ee 
Tie) a+ x | x-+ a | | 2a 
; z => a a ere 5 z | | . | 
Chi gt : 4 | eae wee |sin w,| | : 
TS rb Ae ce ea: | cos? wy, | cos? w, | wes 


Im Fall a) liefern die in Ziff. 3 und 4 angegebenen Formeln folgende Spannungen: 


Fiir x2 < a?: 


Tj = Te = Gs 
(Oy. = (o,)e a elk 
2 2 
()=—(1 + CoP +O toms (36) 
7 2 
(c,),=— 0 cee Co) p — = Cr + O2> 


wobei p und q nach (1’) einzusetzen sind. 
Fir x2 > a?: 
Cito 08 
(o,)1 = (cy). =0, 
2a x 


(o,), = ————— |A sin (Bin 4) + Beos (Binz 4)|-+ Cy + 691> 


us )x2 —a? |s| 


sae ie is 6 
A sin (Bin) + Boos (Bin 24) — Cy + do2- | 


a’ 
— 
wo 
~_ 
—— 


2a a 


Cd = ea Fe 


Die Formeln (36) und (37) wurden in meiner friiheren Arbeit direkt aus den Randspannungen 
p> q berechnet. 
Im Fall b) erhalt man 


c Rf G : 
t= — ee» sin 2 o—E); 


2 2 | 
(oy); + (6,); = ark te eb (n—2 0) (A + Cy tg @) cos @, = c= Cy + 01> (38) | 
(6,)1 — (62) = = (D eh 2) sin 2 @, +. 2 F) —= C— oy, 


wobei 4, B, C, D, E, F aus Tabelle 4 zu entnehmen sind. Die fiir die untere Halbscheibe gelvenden 


Formeln folgen wegen (33) aus (38), indem man f durch — # und w, durch co, ersetzt und das Vor- 
zeichen nach (34) beriicksichtigt. 
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Tabelle 5. Die ldéngs der Geraden x = a (Fall c) geltenden Ausdriicke A, B, C, D, E, F 
fiir die Beanspruchungsarten I bis V 


A 2 bs | ‘ 
Oa — — — I 
1 2 d 2 2 (te o, +6) — 
q : Ms 2 (tg ae | of 
4B Os 1—4p2 Q | < ie Q | Q 
S| Tap = Teen pe a | 2 (tg w, B) = | * 
ao Se ae 
|| ian @) 1+ 4p? @ oe 
Bis 3489 . 0 (1+ 46%) G 
E 
I | Gof 8 (F—«] sin (6 In |sin ol) 
II | Goj £ (3 a 01) cos (f In |sin «]) Bs 
III [ee tg; cos (8 In |sin «,|) sina (5 - on) [4f sin (B In|sin co,|)—(1 42) cos (f In|sin ,|)] Cof Bh = ol ear 
; 2 M/a? 
IV . tg w, cos (f In |sin «,|) Gin B - 0] — [sin (f In|sin @,|) + 2 6 cos (6 In |sin @,|)] Cof B Be o,} i ag 
V p tg w, sin (6 In |sinw,|) Gin B (F _— a1) — [2 sin (f In |sin ~,|) — cos (f In |sin @,|)] Cof B ( = or) (oe 
F 
é P 
74 Sin Bp a on) cos (f In |sin @,|) 
IL | — Sinp (e es or) sin (f In |sin |) © 
Ph ee Qla 
eal —(s tg w,sin(f In|sin w,|) Gof 6 (e ov) -+- [4B cos(f In|sin w,|) + (1—42)sin (6 In |sin @,|)] SinB es o,} 1.4 
| : ’ 2 M/a? 
Vy -P tg, sin (f In [sin @,!) Cof 6 ee oy} [2.8 sin (8 In |sin «,|) — cos (f In |sin @,|)] Ging Me o,} ae 
V 


‘ I 
. tg w, cos (f In |sin @,|) Cof B bog oy) — [sin (8 In |sin @,|) + 26 cos (B In |sin @,|)] Gin B ce o,} G 
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ne | (eek 
[PF 16) s00 ge + (2 © ug) as] fo uy — o(gns)eras— 2nd) svs Ir) 0 12 SOME emda” | 
b+ BIG DS z : 
Paci a nt) wwde— (Eau) so @ gms j ane ( 3) soo ae+(F gu] A glag— 0 ea seeie tI eee tT AI 
pee ) anes dr+tl(oz , 
D/O 3 ug pe vad Di D/O (Serie D Lae t oP eT JP tT 
G nd) gs|%0 dung — (°F a1) ms gs] %0 g og — 0 (08 0 ent 6 avis 
i p (nz ay i ae One te : ; - OS j ; 
ae 5 (Fug) asso gue — (Burg) soot J \ag - aie i 0 II 
aS ) x09 * ; re ‘ GAN fe z BD D 
AEN Ig o gms 2 (°S uf) ws te gop Pare a 0 i: i 
SSS 
a 
fc rc ad a ‘sf aie 5 
a ee 


TH “d “OD “a ‘¥ ayonapsny uapuaya? (p yo) d 


A 819 [ uajapssunyonidsunag arp anf 
SMIPPY Woda sasiasyy uatsojas p = x yung uap wn 


‘uaulayy 4yas sep sSupj aig ‘9 ajoquy, 


S PEP eae IZ—(*m Z—2) oes 0 °0 (cd'% + T) 0 0 9 9Z-| A 
ie eee s | on (oe) leg ee eae nee ‘o $a] 0 ree ort 1 0 |ar 
0 e ; foe Ue 2) J !o9(cd b-1) + (' Z-2) Jing to fr Jy] 0 Pac t—To %) 5 ueeT 0 |zIr 

0 4 2 C0 S— x) J lag 0 58 +o 34) = ON ar 

= (' g—x) ging | 0 ~@ e+ to %)— 0 0 = I 


an SEES eeeemerceeeee ne ee ee ene ee Me a 


A $19 [ uatunssunyonidsuvag a1p anf q “A ‘a ‘0 ‘eK ayonapsnp uapuarjas (q 110.1) asyapr-A sop ssup] ag “py eT[equy, 
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Im Fall c) gilt 


Se cos O, / cos B aes, : 5 : 
Ss a ake < : ~¢ (5 {c sin |f In sin w, + 2 [or = 2 
: 3 GEN) 
=i). Cos t In sin «, 4 5 (o 3) | _ F sin : (0 -| 2| F cos : (or | ap 
4 2 2 2 2) /? 

2 Tae CR, B(=-—o, sels ; Dey, ge) 

(Oy). + (64) = ines, e ( {(4 + 2 C, tg w,) sin p In sin w, + a [os ++ 3) 
¢ 1 > . iE 7 | 2 
+ (B — 2 C, tg w,) cos G In sin w, 4 5 (1 }- zh ts = Cy + 01> , (39) 

Fe oe fos Oy A yc lo tn oi 3 x 
( Or AT Ses ap ae ( D sin |B In sin w, + 5 ( oa 5) 


= G cos lf In sin w, ++ “ (> ae 2)h Sed Syn = (os “fF Z| 


— Ecos s(t 3 ‘\ 


(™ 43) 
A, B, C, D, E, F sind in Tabelle 5 enthalten. 


Fir die untere Halbscheibe braucht man lediglich 6 durch — f sowie w, durch w/, zu ersetzen 
und das durch (34) festgesetzte Vorzeichen zu verwenden. 


SchlieBlich bekommt man im Fall d) 


2 
ee C, aie O01 . 


its 1 SIN M2 2 1, : 3 0 3 
E= eee oes ef @. [Csin (Bing), > 9 02) + D cos (én oat os 0} 
fe | 
Seis Ws 
+ Esin>— F cos FAL 
2 ef : 
(9,). + G1 = — = yz 4 sin G In ae + 2 + Bcos (8 In $5 + || 
a (40) 
2 
= = C, he 001 2 
ea) = B88 Oe een amon’ 2 >. w,\= C cos (fin =o = wr, 
y p es 2 2a 2 2a 2 
a 
: 2 
— sin — Boos | — 2 Gon: 3 | 


A, B, C, D, E, F sind in Tabelle 6 zusammengestellt. Die Spannungen in der unteren Halbscheibe 
folgen sofort aus (40), indem man f durch — 8 und w, durch w, unter Beachtung von (34) ersetzt. 


Die Formeln (40) gestatten die Untersuchung des Spannungszustandes in unmittelbarer Um- 
gebung des Eckpunkts (a, 0). Nahert man sich diesem Punkt auf der Geraden @, = konst., so 
_ sieht man, daB die Spannungen im wesentlichen Terme von der Form 


sin (6 In - + Kons) cos (6 In + kons| 
: und 


yz 2e 
re a 
enthalten. Fiir g —> 0 strebt der Nenner gegen Null; gleichzeitig aber zeigt der Zahler ein eigentiim- 


 jiches Verhalten: Er andert bei Annaherung an den Punkt g = 0 unendlich mal sein Vorzeichen. 
— Dieser Effekt ist allerdings auf eine auBerordentlich kleine Umgebung dieses Punktes beschrankt. 
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| FallI (y=0) 
dias a Pipe Prize 
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a 
Vcr 
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Abb. 6. Spannungsverlauf lings der x-Achse im Fall I fiir C) = 0 und C, = 2/2 (Index I] bzw. 2 fiir die obere bzw. untere Halbscheibe), 


ist * Fall I (w=0) 
———ai _o x BP 
ees H- 0 oe Pm 2a 


Ql 


mae -10 eon a 05 Te 6. 


Fall IE (y~0) 


Gn 2a, 


Abb. 8. Spannungsverlauf langs der x-Achse im Fall ITI fiir C, = 0 und C, = a/2 (Index 1 bzw. 2 fiir die obere bzw. untere Halbscheibe) 
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ieee 2G Faille Tu. (x=0) 
-—- g 


Gn 2a, 


a 


Pees 
— — a ——— 1 — — 
=45. -40 -O5 0 05 40 F rina 


Abb. 9. Spannungsverlauf langs der y-Achse in den Fallen IT und III fiir C, = 0 und Cy = a2. 


Abb. 10. Spannungsverlauf langs der x-Achse im Fall III fiir C, = 0 und C, = m/2 (Index I bzw. 2 fiir die obere bzw. untere Halbscheibe). 


\ 2 (Gcle__, (Om) 
pan 


Abb. 11. Spannungsverlauf langs der x-Achse im Fall IV fiir C) = 0 und C, = a/2 (Index 1 bzw. 2 fiir die obere baw. untere Halbscheibe). 
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Fall IV (x=0) 


Fale 


Abb. 12. Spannungsverlauf lings der y-Achse im Fall IV fir C, = 0 und C, = 27/2. 
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Abb, 13. Spannungsverlauf lings der x-Achse im Fall V fiir C, = 0 und C, = 2/2 (Index I bzw. 2 fiir die obere bzw. untere Halbscheibe). 


: Fall V- (x=0) 


Abb. 14. Spannungsverlauf lings der y-Achse im Fall V fiir C, = 0 und C, = n/2. 
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Betrachtet man beispielsweise die Funktion cos [6 In (9/2 a)], so erkennt man, da8 die hier in Frage 
kommende, am weitesten vom Punkt o = 0 entfernte Nullstelle bei 


ist; fiir den gréBtméglichen B-Wert, namlich 6 = (In 3)/2 2, bekommt man 
=e 3 1,295.10. 


Der oszillierende Spannungsverlauf ist nur vorhanden, wenn fb ~ 0 ist, d. h. wenn beide Halb- 
scheiben aus verschiedenem Material bestehen. Praktisch ist er bedeutungslos, da er in den Bereich 
sehr groBer Spannungen fallt, in dem bei einem wirklichen Werkstoff die Proportionalitatsgrenze 
meist bereits tiberschritten ist, so daB® die lineare Theorie hier nicht mehr zutrifft. 

Um einen Einblick in den Verlauf der Spannungen zu bekommen, wurden diese lings der x-Achse 
und langs der y-Achse numerisch ausgewertet, und zwar fiir die fiinf betrachteten Belastungsarten 
I bis V. Die ausgezogenen Kurven gelten dabei fiir 8 = 0 (d.h. C, = 0), die gestrichelten fir 
B = (In 3)/2 a (d. h. Cy = 77/2); sie stellen die méglichen Extremkurven dar, falls 6 > 0. Ist 6 <0, 
dann vertauschen die beiden Halbscheiben ihre Rolle. Die Spannungen wurden dadurch dimensions- 
los gemacht, da sie im Fall I (Abb. 6, 7) mit dem mittleren Druck p,, = P/2 a, in den Fallen II 
und III (Abb. 8, 9, 10) mit dem mittleren Schub q,, = Q/2 a, im Fall IV (Abb. 11, 12) mit der 
Druckspannung py = (3/2) M/a? (dieser Wert gilt bei linearer Spannungsverteilung bei x = — a) 
und im Fall V (Abb. 13, 14) mit der Vergleichsspannung o, = (2/7) C, dividiert wurden. 


7. Zusammenfassung. Fiir beliebige Punkte eines ebenen Verbundkérpers, der aus zwei ver- 
_ schiedenen, lings einer endlichen Strecke fest miteinander in Verbindung stehenden elastischen 
Halbscheiben zusammengesetzt ist, wird der Spannungszustand vollstandig bestimmt. Es werden 
fiinf Belastungsarten diskutiert, namlich J. reiner zentrischer Druck bzw. Zug, IJ. Schub ohne 
Verdrehungsméglichkeit, III. reiner Schub, IV. reine Verdrehung, V. Warme- und Vorspannun- 
gen. Um einen Einblick in den Verlauf der Spannungen zu bekommen, wurden diese langs der 
Berandung und langs der Symmetriegeraden fiir die beiden méglichen Extremfialle (beide Scheiben 
bestehen aus demselben Material; eine der beiden Scheiben ist starr) und fiir die fiinf betrachteten 
Belastungsarten numerisch ausgewertet und in Diagrammen in dimensionsloser Form niedergelegt. 


(Eingegangen am 27. August 1959.) 


Anschrift des Verfassers: Dr.-Ing. Hans Bufler, Miinchen, Arcisstrafe 21 
Institut fiir Techn. Mechanik der Techn. Hochschule. 
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Zur homogenen Maschine mit Zusatzdrehmassen 


Von H. Bufler und F. Kiefling 


1. Einleitung. Ein Drehschwingungssystem, bei dem samtliche Drehmassen und Federkon- 
stanten einander gleich sind, wird als homogene Maschine bezeichnet. Die Higenfrequenzen und 
Schwingungsformen kénnen in diesem Sonderfall — selbst fiir beliebig viele Massen — explizit 
durch einfache trigonometrische Formeln angegeben werden! oder aus den Nullstellen der Grammel- 
schen Frequenzfunktionen? bestimmt werden. 

Falls an den Enden des homogenen Teiles der Maschine zusatzliche, von den tibrigen ver- 
schiedene Drehmassen und Federkonstanten vorhanden sind, fiihrt die Berechnung der Higen- 
frequenzen bei Verwendung der trigonometrischen Darstellungsform auf transzendente Glei- 
chungen, die numerisch oder graphisch zu lésen sind. Fiir insgesamt N Drehmassen miissen 
sich, wie bekannt, N—1 Higenfrequenzen ergeben. Die Durchrechnung praktischer Beispiele 
zeigt indessen, da® fiir gewisse Zahlendaten dieses Verfahren bei alleiniger Beriicksichtigung der 
reellen Lésungen eine oder mehrere Frequenzen zu wenig liefert. In diesem Fall haben die 
transzendenten Gleichungen fiir die Eigenwerte neben den reellen auch komplexe Loésungen’, 
welche die restlichen — es handelt sich hierbei stets um die héchsten — Frequenzen liefern. 

Das Grammelsche Verfahren ergibt bei Beniitzung der Frequenzfunktionen in der urspriing- 
lichen Form eines Polynoms* des Eigenwertes € simtliche Frequenzen. Wahrend bei der rein 
homogenen Maschine stets 0 < €< 4 ist, kénnen beim Vorhandensein von Zusatzdrehmassen 
durchaus Werte von ¢ auftreten, die gréBer als vier sind. Dieses ist immer dann der Fall, wenn die | 
trigonometrische Darstellung zu komplexen Lésungen fihrt?. 

Da die Grammelschen Frequenzfunktionen nur fiir den Bereich 0 <¢ <4 tabelliert sind, 
andererseits die trigonometrische Darstellung auch gewisse Vorziige besitzt, soll letztere der vor- 
liegenden Arbeit zugrunde gelegt werden. Unter Beschrinkung auf eine oder zwei Zusatzdreh- 
massen werden die Bedingungen untersucht, bei denen komplexe Eigenwerte auftreten, ferner die 
zu ihrer Bestimmung bendétigten Gleichungen sowie die Ausdriicke fiir die Momente und Dreh- 
winkel angegeben. 

Bei den Kurbelwellen von Verbrennungsmotoren hangt die Differenz der Drehwinkel zweier 
benachbarter Drehmassen nicht nur von dem zwischen ihnen wirkenden Drehmoment, sondern 
auch von allen iibrigen Drehmomenten ab. Die Beitrage, welche diese Drehmomente zur Gesamt- 
torsion liefern, werden als Nebentorsionen® bezeichnet. Diese haben auf die Grundfrequenz geringen, 
auf die héheren Frequenzen u. U. wesentlichen EinfluB. 

In nachstehender Untersuchung werden die Nebentorsionen nicht betrachtet, so daB die ge- 
wonnenen Ergebnisse nicht auf Schwingungssysteme angewendet werden diirfen, bei denen diese 
eine Rolle spielen. Mit der trigonometrischen Methode lassen sich aber auch diese Falle erfassen, 
wie in einer demnachst erscheinenden Arbeit der Verfasser gezeigt werden wird. 


2. Grundgleichungen. Bedeuten 0, bzw. g;, das polare Massentragheitsmoment bzw. den Dreh- 
winkel der k-ten Scheibe, c, bzw. M), die Drehfederkonstante bzw. das Torsionsmoment zwischen 


der k-ten und der (k + 1)-ten Scheibe gemaB Abb. 1, so liefert das Momentengleichgewicht an der 
k-ten Scheibe 


M,— M.—1 = O; o (b= Lire IN (1) 


AuBerdem gilt die Formanderungsbezichung 


Hae M, = ¢ Q.41-)): (2) 


P. Funk, Z. angew. Math. Mech. 15 (1935) S. 113. 
Biezeno-Grammel, Technische Dynamik, 2. Auflage, Bd. 2, S. 367ff, Berlin 1953. 
Hierauf geben Th. Péschl und L. Collatz, Z. angew. Math. Mech. 18, (1938) S. 186, einen Hinweis. 
R. Grammel, Ing.-Arch. 2 (1931) S. 228. 
Der Zusammenhang zwischen dem Eigenwert x, der bei der im folgenden verwendeten trigonometrischen 
Darstellungsform auftritt, und dem Eigenwert ¢ bei Grammel ist gegeben durch € = 2 (1 — cos x). 

6 Auf ihre Bedeutung fiir die Berechnung der Eigenfrequenzen von Mehrzylindermaschinen hat erstmals 
R. Grammel, Ing.-Arch, 4 (1933) S. 287 hingewiesen; vgl. auch C. B. Biezeno u. R. Grammel, Technische 
Dynamik, Bd. 2, S. 349, S. 361 und S. 416, 2. Aufl. Berlin/Goéttingen/Heidelberg 1953. 
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Die Elimination des Drehwinkels fiihrt zu 
M, 1 1 M;,_1 My +1 
J -\M, wet WO 
Ck r G | O,, a | ‘ Or; On 44 ie @) 
und mit dem Ansatz 
M, = M,-sinwt (4) 
(@ Kreisfrequenz der Eigenschwingung, t Zeit) 
erhalt man die Rekursionsformel fiir die Momente 
M,,_ 1 1 2\ M,., 
—-—(e +2 =| m, +! =o. (5) 
O;, Op 5 On On 41 


Diese Darstellung, welche auf H. Neuber! zuriickgeht, gestattet gegeniiber der sonst meist ver- 
wendeten entsprechenden Beziehung fiir die Drehwinkel eine einfachere Beriicksichtigung der 
Randbedingungen. Die Rekursionsformel fiir die Drehwinkel, die man aus (1) und (2) nach Eli- 


mination der Momente erhailt, lautet: 


Ce—1 Pk—1 — (& + Ch—1 — w? O,) Y, + Pe 41 = 0- (6) 


M Mm M, My My Myo My. My 


eoeee | \YA-74 (4K fF eocce 


G, O GY 


Abb. 1. Allgemeines Drehschwingungssystem mit N Massen. 


Fiir den homogenen Teil einer Maschine vereinfacht sich (5) wegen 0, = O und c, =c zu 


M,—1—(2— 08 2) My + Mes = 0. (5’) 


Re. 
Diese Differenzengleichung lat bekanntlich die einfache Lésung 
M, = Csin(kx +y):snwt (7) 


(C beliebige Konstante) 
zu, falls man 


2— @22 = 2'cos x, 
also pe 
o S 2 ie Sgn (8) 
iE C) 2 
setzt. Die GréBen x und y sind dabei durch die Randbedingungen festzulegen. Fiir den Winkel- 
ausschlag y, der k-ten Masse gilt 


Pia een 
2esin + 


cos (: =o oo / ‘sinwt, (9) 


wie man durch Einsetzen von (9) und (7) in (2) nachweisen kann. 
Um die Konstante C aus (7) und (9) zu entfernen, ist es zweckmaBig, das Momenten- baw. Dreh- 


winkelverhaltnis M,/M, und 9,/p, zu bilden: 


M;, nee sin (k 4-5 y) (10) 
MM,” sin(x+y) ’ 
cos|(k—)» +7 
ies : (11) 
Pt cos (F + ») 


1 H. Neuber, Ing.-Arch. 22 (1954) S. 258. 
; IN 


ASD) 


H. Bufler u. F. KieBling: Zur homogenen Maschine mit Zusatzdrehmassen Ingenieur-Archiv 


3. Die homogene Maschine mit einer Zusatzdrehmasse. Der homogene Teil der Maschine besitze 
n Massen!, die Zusatzmasse sei mit 9, 1, die dazugehérige Federkonstante mit c, bezeichnet 
(Abb. 2). Fiir den homogenen Teil einschlieBlich der Momente M, und M, gilt (7); die Zusatzmasse 
wird durch (5) mit k =n erfaBt: 
M, 


M e23) | 1 i wo? n+1 
BT sea Se |e Piste M 4 ‘ (ie (12) 
Oo (6 Q;, +1 | can Om. 1 


Aus der Randbedingung M, = 9 (links freies Ende) folgt mit (7) 
ee 
Die zweite Randbedingung M,, ,; = 0 (rechts freies Ende) liefert aus (12), wenn man dort die Mo- 


mente M,, und M,,__; aus (7) einfiihrt und w? durch x entsprechend (8) ausdriickt, die Bestimmungs- 
gleichung fiir die Eigenwerte x 


0) Captian Oeaoe ; — 
(1 | ee A = sin? 4 sinn x —sin(n—l)x=0. (13) 


Abb. 2. Homogene Maschine mit einer Zusatzdrehmasse. 


In gewissen Fallen, d. h. immer dann, wenn im Bereich 0 << « < z weniger als N — | Eigen- 
werte liegen?, existiert neben reellen Lésungen dieser Gleichung auch eine komplexe; mit 


x=a+ip (a,f reell) 
liefert der Realteil der Gleichung (13) den (physikalisch allein sinnvollen) Wert 
C=, 


s 
womit der Imaginarteil von (13) zu folgender Bestimmungsgleichung fir / fiihrt: 


( eee A Soft) Sin mB + Sin (m— fp =O. (14) 
n+1 n 
Die Formeln fiir die Kigenfrequenz, das Torsionsmoment und den Drehwinkel gehen mit y = 0 
und x =a +i iiber in 
wo =2)/£ CofE (15) 
M,, ni1 Cinkp aes 
M, (—1)F+ Sing (ives 1) 25.) (16) 
il 
ein (he 


(b==1,2. 20); (17) 


Der Drehwinkel y, 4.1 der Zusatzdrehmasse errechnet sich aus (6) mit k = n + 1 und mit (15) 


Put _ L 


va 71 nS Pre eae fe 
1—4 6 Gof 5 


(17’) 


* Die Zahl der Massen im homogenen Teil bezeichnen wir immer mit n, wahrend N die Zahl samtlicher 
Massen einer Maschine bedeuten soll. 
* Die Lésungen x = ma (m = 0, + 1, + 2,...) sind trivial, da sie keine Schwingung liefern; auBerhalb 


des angegebenen Bereichs liegende Eigenwerte fiihren wegen der Periodizitat der Gl. (13) zu den schon vorher 
erhaltenen Eigenfrequenzen. 
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Zwecks Diskussion der Gleichung (14) formen wir sie folgendermaBen um: 


B C, (Sin (n —1)B 6) 
Pyle Ns 
Co} BLE ( Sinnp ie 6) : (18) 
SRS heats TEE n+1 
L R 


Die durch jede Gleichungsseite gegebenen Funktionen L — L(p) und R = R(f) lassen folgende 
Schliisse zu: 


1. Ist die durch die Ungleichung 


il 
a (19) 


gegebene Bedingung erfiillt, so liefert (13) N— 2 reelle Eigenwerte und (14) oder (18) den rest- 
lichen (héchsten) Eigenwert!. 


2. Ist (19) nicht erfiillt, so liefert (13) samtliche N — 1 Eigenwerte [Gleichung (14) hat in diesem 
Fall keine reelle Lésung fir £]. 


4. Die homogene Maschine mit zwei Zusatzdrehmassen auf einer Seite. Der homogene Teil der 
Maschine bestehe wieder aus n Massen; am rechten Ende kommen die Zusatzmassen O05 
und @,, 2 und die Federkonstanten ¢c, und ¢, . 1 hinzu (Abb. 3). Fiir die Momente im homogenen 


M M, Mp Mn-1 My Mye7 My+e 


G, +2 


Abb. 3, Homogene Maschine mit zwei Zusatzdrehmassen auf einer Seite. 


Gre 


Teil (einschl. M, und M,) ist (7) zustandig; die beiden Zusatzmassen werden durch die Rekursions- 
formel (5), angeschrieben fir k =n und k =n + 1, beriicksichtigt: 


M1 1 1 wo M41 
can at Mee | M. ft) 9 20 
iC) (6 w On = i a OO 3 ( ) 
M. 1 1 wo? M,, 4.2 
areas | Ni pee ee (21) 
ar n+1 
Ona bean O42 Cnt ;) O,42 


Die Randbedingung am linken Ende M, = 0 ergibt mit (7) 
re Wy, 


Die Randbedingung am rechten Ende lautet M,,2 = 0. Eliminiert man nun M,+1 in (20) 
mittels (21) und verwendet fiir M,,_; und M,, die Ausdriicke (7), so folgt mit (8) 


@) oe O Cree 2a% (0) 2) c ck 2) 
sey nee ay 2 2 —4 sin? -— || sinn x 
(6 ( 2 OTe Cn ron aloe + O,42 Cea 2 


n+1 
(0) (0) Cs OAS aN 
oo RE a ae sin? )sin (x —1) 2 =0. (22) 
oe ar On +42 Cn +1 2 


Zur Exmittlung von etwa vorhandenen komplexen Lésungen dieser Gleichung setzen wir, wie 


in Ziff. 3, 


x=a+ip (o,f reell) 
und erhalten 
ot 


1 Wegen der Symmetrie der Kurven L = L(f) und R = R() beziiglich der Achse 8 = 0 ergeben sich Zwei 
Werte fiir B die ck nur durch das Vorzeichen unterscheiden, jedoch beide zum selben Ergebnis fiihren. 
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Sin n fp 


fatale ai a oka ta 


n+1 n+2 Creel 


-(g-+e5-t*s CoP) Sin (nm — YA =0- C3) 


n+1 n+1 


n+2 


Fiir w sowie fiir M, und y, im homogenen Teil gelten hier die Beziehungen (15), (16) und (17). 
Das Moment M, 1 kann entweder aus (20) oder (21) berechnet werden; aus (21) ergibt sich 


oe 
ea 
Meany te aye 7 n+1 5 Sinn fp ; (16’) 
Wee SAO ee he} Ac oor 2 Sin B 
Oo a Oras Cn +i = 


Fiir die Drehwinkel gn 41 und @n+2 der Zusatzdrehmassen bekommt man aus (6) mit (15) 


De B 
yee fon a Gof? > 
Pn+1 Pn : €n 1 () : 2 (ia 


: a 5 (o) 
Pi Py Ca 1 oO: ee tC, Pe B ( Cc n+2 72 B 
Sof? 1—4 1) 1 
f Cn Gioe Co} 2, Cn 1 O 2 


Pn +2 =) Pn +t ute 1 at eg (17’”) 


z & 
V1 Pi ey oD n+2 Gof? B 


Cn+1 C) 2 


Zur Diskussion von Gleichung (23) ist folgende Form giinstig: 


(0) 2 
eeaB (o,, : Sin (n — 1)8 C) Sy 
SURG Se aes Oo ee ak = ete PL ae 2m 
Ory @, 12 CL, £1 Be iy 


L R 
Die Betrachtung der durch jede Gleichungsseite gegebenen Funktionen L = L(f) und R = R() 
fiihrt zu den Ungleichungen 


) 2) c : 
ees ee 29 
O55 a O, 42 en ( ) 
() 2 
ce () pao 1 

is 3 26 
Aa o... 6 O eee me i e) 

Sera Or .3 Cree 


und es gilt: 
1. Falls die Bedingungen (25) und (26) erfiillt sind, liefert (22) N—3 Eigenwerte und (24) 
die restlichen zwei (es handelt sich hier um die beiden héchsten Eigenwerte)!. 


2. Falls eine der Bedingungen (25) oder (26) erfiillt ist, liefert (22) N — 2 Eigenwerte und (24) 
den restlichen héchsten Eigenwert. — 


3. Falls weder (25) noch (26) erfiillt ist, liefert (22) samtliche Eigenwerte [Gleichung (24) hat 
in diesem Fall keine reelle Lésung fiir £]. 


5. Die homogene Maschine mit je einer Zusatzdrehmasse auf heiden Seiten. Die Zahl der im 
homogenen Teil der Maschine gelegenen Massen sei wieder n; die linke bzw. rechte Zusatzmasse 
und Federkonstante seien mit @, und c, (bzw. On 41 und c,) bezeichnet (Abb. 4). Fiir die Momente, 


1 Sinngema8 gilt auch hier FuGnote 1 von Seite 253. 
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die im homogenen Teil liegen (einschl. M, und M,), gelten die Beziehungen (7), fiir die beiden Zu- 
Satzmassen die Formel (5) (und zwar mit k = 0 und k = n) 


M_, 1 tke M, 
Op ac | o— a) Mo agp (27) 
M 
n—l1 1 1 Ww? M 
— — | oa 8 n+I1 
a) (6 On41 =| Ma iy omer (28) 


On 41 


Abb. 4. Homogene Maschine mit je einer Zusatzdrehmasse rechts und links. 


Damit und mit M,, M,, M,—1 und M,, aus (7) sowie w? aus (8) erhalt man 


: C) Char ek : 
sin x cos y + eos & (1 6, 4 Fs sin? >) sin y= 0, (29) 
sin (n— 1) x —{1 4 as 4 © sin? * | sin n x| cos y 
0, +1 Cn 2 } 
2) Cea X i 
-|- cos (n—1)x (1 == eine | cos n 7 sin'y = 0. (30) 
\ On on 2 


Eliminiert man daraus y (was gleichbedeutend ist mit dem Verschwinden der Koeffizientendeter- 
minante vorstehender linearer homogener Gleichungen fiir sin y und cosy), so gelangt man zu 
folgender Bestimmungsgleichung fiir die Eigenwerte: 


t 7" 4 sin? alt +o 4G sint 5 | sin mx 


0 n+1 n 


[it Seg) (t+ 8 
n+ 


4° sin? _ sin (2 1) 2 sin (n — 2) — 04. (31) 


Zur Untersuchung der komplexen Liésungen dieser Gleichung setzen wir wieder 


=~ 1p (x, B reell) 
und bekommen 
70 
sowie 
(0) c ; 
fi gst) (of, 18) ane 


n+1 


+. (0 [se ays 2 af ( epee ey Gai zen (n—1) B + Gin (n—2) B =0. (32) 
(OF Co 2 Oe G, 2 
Die GroBe y berechnet man am einfachsten aus (29); mit x =a +1 folgt 


Gof B + ( + 2 4 = oop 8) 


ctgy =ibigiy =1 Sink , 


also ein rein imagindrer Wert fiir y. 
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Fiir die Kreisfrequenz gilt wieder (15); die Momente und Drehwinkel im homogenen Teil 
werden aus (10) und (11) mit x =a +i und y aus (33) ermittelt: 


) c Bese 
Sin (k — 1 -{1- 4 — Wp}? — ein p 
wale 9 pa Sin ( . De (14 Og eEees 2) ez (k = 0, L 9. any (a) 
Mee Se SF cose C) ; 
at ( tet om 2) Sing 
: ; 3) ( C) ees 4 Res Pipe 
Sin | k —- - {1-4 4 — Sof? —] Sin |k p 
Re clos 52 ee a Se 2S oe eee 
C (C) c Neue 
? yt en Pee 
1 (3; 4 e oj 2] ins 
Fiir die Drehwinkel g) und ¢, 41 der Zusatzdrehmassen folgt aus (6) bei Beachtung von 15) 
PO ee oe (35’) 
Yi c O Fpi2 D 
Lane & @O P} 2 
und 
Pa+h Pn = = Eafe ME? ue (35’’) 
al Pi c n+l 2 Bp 
: ae 6) Gof 2 


Gleichung (32) liefert héchstens zwei reelle Werte fiir /. 
Die Diskussion dieser Gleichung sei auf den vereinfachten Fall gleicher Zusatzmassen @,, 41 
= @, und gleicher Zusatzfederkonstanten c, = cy beschrankt. Dafiir gilt 


( +9, — 42605 ) Gin np ! 2(1 Le 4 = Go? 5 | Gin (n — 1) B + Gin (n— 2) B =O. 


Lést man diese Gleichung nach (1 oo se —A — Gof? a auf, so erhalt man nach kurzer Umformung 
die iibersichtlichere Form : 3 


YC eee © , Sinm—p , Ging o 
Seber gee “A 6, = <Cinn heals (36) 
—S=<S — 

L R 


Die Untersuchung der durch jede Gleichungsseite gegebenen Funktionen L = L($) und R = R(§) 
zeigt: 
1. Ist die Bedingung 
[Wek apd ype (37) 


erfiillt, so liefert (31) N—3 Eigenwerte und (32) die restlichen beiden (héchsten) Eigenwertel. 
2. Gilt die Bedingung 

c= O 

9 Oo 


2— Lc 4 <2, (38) 
so erhalt man aus (31) NN — 2 Higenwerte und aus (32) den restlichen (héchsten) Eigenwert. 
3. Ist 


4———_ > 2, (39) | 
so liefert (31) samtliche Eigenwerte; (32) hat in diesem Fall keine reelle Lisung fiir f. 


6. Beispiele. Wir legen den drei Beispielen eine Sechszylindermaschine (n = 6) zugrunde. 
Im Fall einer Zusatzmasse (Beispiel 1) wahlen wir c/ce, = 1/3 und 0/0,=1/10. Durch Ein- 
setzen in (19) sieht man sofort, daB diese Bedingung erfiillt ist. Der héchste Eigenwert muB in- | 


folgedessen aus (14) berechnet werden (vgl. hierzu die graphische Darstellung Abb. 5); Gleichung (13) 
liefert namlich nur fiinf Frequenzen. ‘ 


1 Sinngema® gilt auch hier FuBnote 1 von Seite 253. 
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Wir erhalten 
B = CON 


und damit aus (15), (16), (17) und (17’) 


Ghee 5156 Hee : 


M;, k+1 Gin 0,791 k 
“Shae Sy eae | SMO ie Pe __ k +1 Gin [0,791 (k —0,5)] 
M, co 0,876” 1 or 0,406 tee: 
#2 = 4,420), 
Pi 


Fir die fiinf ersten Eigenwerte ist (13) maBgebend. Auf eine diesbeziigliche zahlenmafige Aus- 
wertung wird verzichtet, da diese wie iiblich geschieht. 


-42 “U8 —O# 0 OY 08 42 


Abb, 5. Ermittlung des héchsten Eigenwertes fiir das Beispiel 1. 


Cicicicici{c, C 
20 z 


000000 8 2B, 


=94 0 4 08 12 16 


Abb, 6. Ermittlung des héchsten Eigenwertes fiir das Beispiel 2. 


Im Fall von zwei Zusatzmassen auf einer Seite (Beispiel 2) legen wir die Daten 
c/eg = 1/3, c/eg = 1, O/O, = 1/10, O/O, = 1/50 zugrunde, welche die Bedingung (26) erfiillen, 
(25) dagegen nicht. Nach dem unter Ziff. 4 Gesagten liefert (22) nur die ersten sechs Kigenwerte, 
wahrend der siebente nach (23) oder (24) ermittelt werden muB. Abb. 6 zeigt die graphische Be- 
stimmung von / entsprechend (24). Es ergibt sich 


[8] = 0,793 
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und damit fiir die héchste Eigenfrequenz und die dazu gehérigen Momente und Drehwinkel 


wo, = 2,159 /< , 


M, _ Die Sin 0,793 k Pe __ (— i +1 Gin 0,793 (k — 0,5) (hs 199 Soy 
Mies ey OB 19 see Py 0,407 

a ees or Oro G8 2 5 6.206% 

Maes 1,461 , = Sone oe 


Im Falle von je einer Zusatzmasse rechts und links (Beispiel 3) wahlen wit ety = c/Cg 
= 1/3, 0/0, = O/O, = 1/10. Mit diesen Daten ist die Bedingung (37) erfullt. Nach Ziff. 5 liefert 


dann (31) die ersten fiinf Eigenwerte, waihrend die restlichen beiden (auf die wir uns hier beschran- 


eee 
ColCICICICIC|Ce|| % & 3 


4 0, 900000 O, 


16 “48 ~G8 ~O4 0 OY 08 12 46 


Abb. 7. Ermittlung der beiden héchsten Eigenwerte fiir das Beispiel 3. 


ken) nach (32) bestimmt werden miissen. Die entsprechende graphische Darstellung gibt Abb. 7. 
Man entnimmt daraus 
0,755 
B| = 


0,805 


und erkennt ferner, da die Differenz dieser beiden Werte um so kleiner wird, je gréBer diese sind. 
Fur die entsprechenden Frequenzen bekommt man? 


W, = 2,152 ee : (_ = 2,164 as 
Die Gleichungen fiir die Momente und Drehwinkel lauten hier nach (34), (35), (35’) und (35”) 
Gin 0,775 (k — 1) — 0,444 Gin 0,775 k 


} : TEMS 
Mere eats pee eee ae 
M, Gin 0,805 (k — 1) — 0,461 Gin 0,805 k fi 
Si 0,413 eee 
Sin 0,775 (k — 1,5) — 0,444 Gin 0,775 (k — 0,5) “ 
~ it CO fir, 
Oe 4 4) : (k ==T; 2). . 6), 
1 Sin 0,805 (& — 1,5) — 0,461 Sin 0,805 (k—0,5) 
0,604 Ur @, - 
eles 0,0693 fiir m,, 7 — 0,0693 fiir we, 
%:  \—0;0685. far a, % | 0,0685 far w,. 


+ Das von Th. Poschl und L. Collatz a. a.O. fiir eine Sechszylindermaschine mit zwei gleichen Zusatz- 
massen und gleichen Zusatzfederkonstanten angegebene Nomogramm fiihrt mit unseren Zahlendaten zu 
Ws = 1,90 Ve/ 0, w, = 1,96 Ve/ ©. Anscheinend sind in diesem Nomogramm die Parameter y und y’ ver- 
wechselt, denn nach Vertauschung der beiden erhalt man Ubereinstimmung mit unseren Ergebnissen. 
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7. Zusammenfassung. Unter gewissen Bedingungen liefern die fiir die Berechnung der Eigen- 
werte einer homogenen Maschine mit Zusatzdrehmassen maBgebenden transzendenten Gleichungen 
bei alleiniger Beriicksichtigung der reellen Lésungen zu wenig Eigenfrequenzen. Die fehlenden 
Frequenzen — es handelt sich hierbei stets um die héchsten — erhalt man aus den komplexen 
Lésungen. Es werden, unter Beschrankung auf eine oder zwei Zusatzmassen, die Bedingungen 
untersucht, bei denen komplexe Eigenwerte auftreten, und die in diesen Fallen geltenden Gleichun- 
gen fiir die Kigenfrequenzen, Momente und Drehwinkel abgeleitet. Neben dem Verhaltnis der 
Tragheitsmomente ist auch dasjenige der Federkonstanten fiir das Auftreten komplexer Lésungen 
von Bedeutung, so da®B solche, wie man aus den zustandigen Bedingungsgleichungen ersieht, 
auch bei sehr grofen Zusatzmassen vorkommen kénnen. Man kann iibrigens zeigen, daB beim 
Vorhandensein von z Zusatzmassen maximal z komplexe Eigenwerte méglich sind. 

Es sei hier noch vermerkt, da sich die gewonnenen Ergebnisse analog auf homogene Schwin- 
gungsketten mit Zusatzmassen an den Enden itbertragen lassen. 


(Eingegangen am 7. Oktober 1959.) 


Anschrift der Verfasser: Dr.-Ing. H. Bufler und Dipl.-Ing. F. Kiefling, Miinchen, Techn. Hochschule, 
Institut fur Technische Mechanik. 
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On the steady thermoelastic problems with elastic characteristics 
depending on temperature 


By M. M. Stanisic and R. M. McKinley 


1. Object of Investigation. A comparison is in order between the solution of the steady-state, 
axially symmetric stress field presented by Trostel! (hereafter referred to as an exact solution A) 
and the approximate solution presented by Stanisié and McKinley? (hereafter referred to as 
solution B). Trostei determines a solution in series for the Navier equations and employs a two- 
term truncation for actual calculations. In solution B, Stanisié and McKinley replace the com- 
patibility equation, in terms of a stress function, by an approximation and solve this by quadra- 
tures. The object of this paper is to show that the technique used in solution B is equivalent to 
a two-term truncated solution in series of the compatibility equation. The similarity between 
solutions A and B will elucidated in this study. 


2. Analogy between two Solutions. To begin the analogy, consider first the case for which. 


Pi: =Po= (no external load) , 
o(T) =a[T(r)]. 
ET) ==, (constant) . 
Equations (18), (19), and (20) of solution B are evaluated as 
50) jp 
(a) ®,(r) = 5 r2 ? 
T(r) a; ee T(é) 
(b) br) = B, | { (0) 40 —(7)' | (0) 40 — = | EI. J (0) 40) a8}, 
0 ) r 0 
nat ac ) 
(c) By) =, E a | 2 | a(@) dO | +4 | é | x(@) dO a 
0 Ts 0 


The stress field is given by equations (29) and (30) of solution B 


= Ey P3(ro) | Pa(r) @,(r) 
(d) CS t=» Eo D,(r)) 3.) ry 
E(T) 
(e) ia Ey . ®3(ro) ®,(r) s ®,(r) a anes he + D,(r) 
; eet P,(r9) P,(ro) 


Equations (d) and (e) are identical with equations (30a) and (30b) of Trostel’s solution A, and both 
exact solutions. 


For the similar case, except E( T) is to be regarded as a variable, we integrate equations (18) 
and (20) by parts to obtain 


2 2 2 
Lei 


(a) D0) = yo ae tee | eee 


1 R. Trostel, Ing.-Arch. 26 (1958) p. 416 (Dezember, 1958). 
> M. M. Stanisié and R. M. McKinley, Ing.-Arch. 27 (1959) p. 227. 
3 M. M. Stanisié, Lectures in Thermoelasticity, Division of Engineering Sciences, Purdue University, 
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Re = T(&) r Te 
3 £ og 
(b’) @,(r) = 7 aie | o(@) d@| dé —2¢ i 8) dé +9 ET)" | 2(@) dO 
g if 
r 0 Le 0 


r é T(n) 


=p [ °F les | «40 dn\ dé 


Pe oe. #7 1 T@) 
= 2e | DE) es | t | 0(@) dO | dé \dy\ dé 
ine Tie Ts 0 


u Vv v as / 
where* 


eee L) 


Collect terms with and without ¢-coefficients, and, also, in equations (d) and (e) of this paper 


replace E, by the variable E(T), denoting the results by 0,9 and dg5. Then the stress field of 
_ solution B can be written in the form 


(f) O, = 07,0 +EOr1, 
(g) 00 = 90,0 + £061. 


This establishes the analogy between solutions A and B. In equations (f) and (g), the o,,9 and 
Go,0 are identical with Trostel’s “Grundlésung” as given by equations (30a) and (30b) of A. How- 
ever, our first corrections o7; and 06,1 are not the same as the analogous o,,; and 06,1 deduced 
by Trostel in equations (35a) and (35b). For example: 


OF) = OF 1— & 5 
where 
Pee (O60) aT 22 PL Fo aT | 
ae ¢ SULA ey tlle 1) 
Se ss, fe ie E(n) dy 01.0) dy Obes Paes pf: | Ge dy o,.0)dn der 


~ 


that is, our first corrections differ from Trostel’s in that they neglect the term in his equations 


(35a) and (35b) with coefficient ae . This term constitutes less than 30% of the total corrections, 


3. Conclusion. We therefore conclude that the approximate solution B gives in a single operation, 
the “Grundlésung” and about 70% of the first correction term of the exact solution A. 


(Eingegangen am 20. Oktober 1959.) 


Anschrift der Verfasser: Prof. Dr.-Ing., Dr. Phil. Milomir M. Stanisié, Division of Engineering Sciences, 
Purdue University, West Lafayette, Indiana; 


Graduate Student R. M. McKinley, Chemical Engineering, Purdue University, West Lafayette, Indiana. 


* The ¢ is a constant. This notation is in agreement with that employed by Trostel. 
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Knickung verwundener Stabe unter Druck einer konservativen, 
kontinuierlich und gleichmafig verteilten Belastung 


Von H. Leipholz 


1. Problemstellung. Als Voraussetzungen sollen gelten: 


1. Die Stabe seien gerade, elastisch, lang und diinn, und sie mégen unter der Belastung nur 
kleine Verformungen erleiden. Dann gelten fiir Gleichgewichtsbetrachtungen am herausgeschnit- 
tenen Stabelement die Gleichungen von Kirchhoff und Clebsch, und man kann iiber die Verkniipfung 
der resultierenden Spannungsmomente in den Stabquerschnitten mit der Kriimmung und Ver- 
drillung des Stabes einfache Aussagen machen, wie sie im folgenden benutzt werden sollent. 


2. Die Stabe mégen fiir einen Querschnitt zwei verschiedene Haupttragheitsmomente beziiglich 
des Querschnitts-Schwerpunktes haben, die aber iiber die Lange des Stabes konstant seien. 


Ferner sei der Stab verwunden, das heift, die Schwerpunkts-Hauptachsen seien schon im 
spannungslosen Zustand um einen konstanten Winkel t je Langeneinheit des Stabes gegeneinander 
verdreht. 


3. Die Belastung der Stabe sei durch kontinuierlich iiber ihre Achse verteilte Druckkrafte q 
gegeben. Dabei mége die Wirkungslinie von q richtungstreu bleiben, das heiBt, stets die Richtung 
der unverformten, geraden Stabachse beibehalten. Dann ist die Belastung konservativ, und es 
kann zur Ermittlung der Eigenwerte des Problems das statische Stabilitatskriterium und damit 
die Gleichgewichtsmethode benutzt werden?. 


Diese Arbeit befaBt sich mit der Weiterfiihrung eines Problems, das von E. Hui? fiir den ver- 
wundenen Knickstab mit Einzellasten an den Stabenden gelést wurde, und es steht auBerdem in 
nahem Zusammenhang mit Arbeiten von H. Ziegler’ iiber die gedriickte und tordierte Welle. 


Mit Riicksicht auf die Lagerung der Stabenden werden die fiinf Knickfalle von Abb. 1 unter- 
schieden. 


2. Aufstellung der Differentialgleichungen. Unter Anwendung der Gleichgewichtsmethode erhalt 
man durch Betrachtung eines herausgeschnittenen Stabelementes (Abb. 2) die Vektorgleichungen 


dg ; dy 5 
eo Aa 4a, tix K=0, (1) 


wobei & den Kraftvektor, Nt den Momentvektor, { den Vektor der Belastung, t den Einheitsvektor 
in der Tangentenrichtung der Achse und s als unabhangige Variable die laufende Lange der Stab- 
achse bedeuten. 


1 A, E. H. Love, — A. Timpe, Lehrbuch der Elastizitat, S. 446, 447. Leipzi i 
; ; i , 5. 446, ; pzig und Berlin 1907. 
> H. Ziegler, On the Concept of Elastic Stability, aus: Adva in Applied M i 
SE Hi Oster Tg Aceh, OS) Coon y vances in Applied Mechanics, IV, New York 1956. 
* H. Ziegler, Z. angew. Math. Physik 2 (1951) S. 265 und 3 (1952) S. 96. 
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Diese Vektorgleichungen seien auf das Hauptachsensystem des Stabes bezogen (Abb. 3). Dieses 
System ist beweglich. Sein Ursprung fallt jeweils in den Schwerpunkt des gerade betrachteten 
Stabquerschnittes, gleitet also entlang der Stabachse. Die €- und n-Achsen decken sich mit den 
Schwerpunkts-Hauptachsen des jeweiligen Stabquerschnittes. Die ¢-Achse fallt mit der Tangente 
an die elastische Linie im betreffenden Querschnitts-Schwerpunkt zusammen. Die Achsen sind 
so orientiert, daB man ein Rechtssystem (€, 7, ¢) hat und die ¢-Achse im Sinne wachsender Lange 
der Stabachse zeigt. 


Mam 


Abb. 2. Abb. 3. 


Da der Stab verwunden ist, fiihrt das Hauptachsensystem auch noch eine Drehung aus, selbst 
wenn man von der geringen zusatzlichen Verwindung infolge elastischer Deformation absieht. 
Der Vektor tv dieser Drehung je Langeneinheit der Stabachse ist mit der elastischen Deformation u 
des Stabes verkniipft. Diese habe beziiglich des Hauptachsensystems die Komponenten 


a) 
v= (2) 3 (2) 
Me 


Dabei sind wz und @, die Kriimmung und 7; die Verwindung des deformierten Stabes. Unter der 
getroffenen Voraussetzung kleiner Verformungen soll insbesondere tT; = 0 gesetzt und die genannte 
Verkniipfung gegeben werden durch - 


0 ‘ We 
os (°) tue (o) ; (3) 
T US 


Schreibt man das Gleichungssystem (1) vom Standpunkt eines mit dem bewegten Hauptachsen- 
system fest verbundenen Beobachters, so erhalt man 


d& : dyn b-@aa 
G, + wxt=—E, ae t+ wxM+ix =o. (4) 


In Komponenten beziiglich des Hauptachsensystems geschrieben gibt (4) die Gleichungen von 
Kirchhoff und Clebsch, hier allerdings in etwas erweiterter Form wegen des Auftretens der kon- 
‘stanten, gleichmaBig verteilten Belastung f: 


dM a 
408 +N; — Qt + Pe = 0; = + T;o,—M,t— Q, = 0, 
dM, a 
oer + Ost — N; 0% + Py = 9, Speak le rae On Ys (5) 
dT; as, 
ONE 4-Q, 0p —Qp Oy + 4 ae, Beat Sr: 


Dabei sind folgende Bezeichnungen fiir die Komponenten gewahlt: 


. M. 0 (Pé 
Nr T; a q 
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Sind a und b die konstanten Biegesteifigkeiten, sowie c die konstante Torsionssteifigkeit des 
Stabes, so ist nach Voraussetzung 1 von Ziff. 1 der Zusammenhang zwischen )¢ und u 


Z M; =a, M, =6ao,, T,=ct =0. (7) 


Die Gleichungen (5) und (7) wiirden unter Beriicksichti- 
gung entsprechender Randbedingungen zur Lésung des 
Problems geniigen, wofern die Komponenten pz, p, und q- 
der Belastung f des Stabes beziiglich des Hauptachsen- 
systems bekannt waren. Gerade das wird aber iblicher- 
weise nicht der Fall sein, sondern es wird die Stabbelastung 
wohl meist mit Komponenten beziiglich eines raumfesten 
kartesischen Systems gegeben sein. Es ist daher noch zu 
zeigen, wie man Vektoren vom Hauptachsensystem auf ein 
rauinfestes System umrechnen kann. 

Abgesehen von der belanglosen, gleichformigen Transla- 
tion entlang der Stabachse mit der Geschwindigkeit 1 
(wofern man die laufende Koordinate s der Stabachse 
jetzt als Zeit deutet) fiihrt das Hauptachsensystem noch 
eine Kreiselung mit 0 um seinen Ursprung aus, den man 
mit dem Ursprung des raumfesten Systems zusammenfallen 

pe lassen kann (Abb. 4). Fiir die Umrechnung von einem 

System ins andere kann man daher Beziehungen verwenden, 

die in der Kinematik des Kreisels gelaufig sind. Unter Verwendung der Eulerschen Winkel erhalt 

man fiir die Umrechnung eines Vektors a von seiner Darstellung a\#45) im Hauptachsensystem 
zu seiner Darstellung a(’f) im raumfesten System die Transformationsformel 


q(H#AS) — Of arf), (8) 


Dabei ist Yf die folgende Matrix: 
cos y cos p — sin y sin @ cos # sin y cos@ + cosy singcos# sing sin? 
Y =| — cos y sin p — sin y cos pm cos} —siny sing + cosycospcos? cospsin# |. (9) 


sin y sin } — cos y sin # cos # 
Ferner gelten fir v noch die Beziehungen 
Oz =cosph+sing sind», 
oO, =—sing d+ cos sind p, (10) 
t=costy+o@. 
Punkte sollen dabei die Ableitungen nach s bezeichnen. 


Wegen der Voraussetzung | kleiner Verformungen kann # nur ein sehr kleiner Winkel sein und 
man kann cos #? & 1 setzen. Dann erhalt man 


tTRey+p, alo pto=y=f[rds=rts (11) 
0 
und damit statt 2{ 
/ cost s sinT s sin ~ sin ? 
Ate cost s cos y sin # }. - (12) 
siny sin? —cospsin ? 1 


_ Fihrt man nun den lotrechten Einheitsvektor ¢ ein mit den Komponentendarstellungen 


Cs 0 
e(HAS) = ey bzw. elf) = (0) A 
€¢. 1 


so erhalt man aus (8), wenn man noch statt 2 das Q[, verwendet, 


e¢ = sing sind, e, = cosp sin}, e=l1. (13) 
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Mit (13) erhalt man statt Q{,, wenn man fiir die dritte Zeile auBerdem sin 0 ~ 0 setzt, 


costs sints ¢ 
YU, =|—sints costs e,). (14) 


Ose 0 ] 


Mit dieser Matrix Q, soll kiinftig die Umrechnung von einer Darstellung eines Vektors in seine 
andere Darstellung gemaB 


qty ==, avs) (15) 
erfolgen. 


Die Gleichungen (13) gestatten noch, zusammen mit der Voraussetzung cos } ~ 1, eine Um- 


rechnung der Gleichungen (10). 


Man multipliziere die erste Gleichung (10) mit cos 9 ~ 1, addiere und subtrahiere gleichzeitig 
rechts den Ausdruck ¢ sing sin 9; dann kommt 


Ws =— sing sin? + cospcos 0) + sing sin 0 (cos 0p +). 
Aus (13) hat man aber 


ee = sing sind, 


e, = — sing sin & + cos y cos od 
und aus (10) noch 
tT=cos¥y+@. 
Also wird 
We =e, + T ep ) 
und entsprechend (16) 
WD; = — Cs Tey 


Die weitere Rechnung erfolgt unter Verwendung von (5), (7), (15) und (16). 


In einem raumfesten, kartesischen System (x, y, z) habe St die Komponenten 


Q. 
RK =] Q, 
N, 


z. 


8 


Den ersten drei Gleichungen von (5) entspricht im raumfesten System die Vektorgleichung 


und f, wie vorausgesetzt, die Komponenten 


d&® 


Se eee is 
oder in Komponenten 
dQ, dQy dN ia oe 17) 
i re ae eae Ca ( 


Das System (17) 14Bt sich sofort integrieren. Seine Lésung lautet 
At 
cs eased & (18) 


wobei A* und B* Integrationskonstante sind. < 
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Man hat damit eine Darstellung von § im raumfesten System. Die Darstellung von % im 
Hauptachsensystem erhalt man daraus gema8 (15) zu . 
* 
Q: costs ° sin7T Ss es 


B* 
RHAS) —| Q, ]=|—sints costs e, 3 
N; 0 0 (3 —fqds 
also 
Q: / A* costs + B* sint s—e q(l—s) 
Q, |= — A* sints + B* cost s —e, q(l— 5) }. (19) 
Ne Tie) 


Jetzt werden die drei letzten Gleichungen von (5) behandelt. Zuerst wird (7) eingefiihrt und 
gibt 


dwe 


a —-— brow, —Q, =9; 
poet | ate t+ Or =O; 


ds 
b Ww, Wg —40,0,=90. 
Wegen der Voraussetzung kleiner Verformungen sind @¢, w, kleine GroBen, und das Produkt w; w, 


kann daher als GréBe, die von héherer Ordnung klein ist, vernachlassigt werden. Es verbleibt somit 
ein System von nur zwei Gleichungen, in das jetzt (19) eingefiihrt wird. Das gibt 


af brew, + A*sints— B* costs +e, g(l—s) =0, 
boon + at wg + A* costs + Bt sin Ts — eg g(JI— s) = 0. 


Verwendet man auch noch (16) und ordnet um, so hat man 
ae, + (a + b)t eg + [g(}— s) —b7*] e, = — A* sint s + B* Bae 


20 
be,— (a + b) te, + [q(J— s) — at] eg = A* costs + B*sints. am) 


Das vorgelegte Problem ist damit auf die beiden Differentialgleichungen (20) zuriickgefiihrt. 


3. Die Randbedingungen. Die Randbedingungen des Problems ergeben sich aus der Lagerung 
der Stabenden. Zur Aufstellung der Randbedingungen sind daher die Knickfalle gemaB Abb. 1 zu 
betrachten. Dort kommen folgende Lagerfalle vor: 


a) Die Stabenden sind frei verschieblich und damit querkraftfrei. 

b) Die Stabenden sind gelenkig gelagert und damit frei von Biegemomenten. 
c) Die Stabenden sind fest eingespannt. 

d) Die Stabenden sind unverschieblich. 


Aus a) folgt 
Q.=A*=0, 0, = Bt =0, (21) 


so das die Gleichungen (20) homogen werden. 
Aus b) folgt fiir s = 0 oder s = | 


M; =ao,=0, was nach (16) gibt e, tre = 0 : (22) 


M,=bo,=0, was nach (16) gibt —e,+te,=0. 
Aus c) folgt: an der Kinspannung fallt der Tangentenvektor der Stabachse t mit dem lotrechten 
Einheitsvektor ¢ zusammen. Im Hauptachsensystem hat t die Komponentendarstellung 
0 
t={0}], 
1 


ee 


XXIX. Band 1960 : H. Leipholz: Knickung verwundener Stibe 267 


womit sich, weil e = tf gelten soll, fiir das eingespannte Stabende im Hauptachsensystem ergibt 


eg 0 
e=j]e,|}/=|0 
e i 


Die Randbedingung heift also in diesem Fall fiir s = 0 oder s = 1 


ee—0, e, = 0. (23) 


Aus d) folgt: der vom Stabanfang zum festgehaltenen Stabende gezogene Ortsvektor muB die 
Horizontalprojektion > 
v = = 
ad AY 


wobei t, ein Vektor mit den Komponenten x und y beziiglich des raumfesten Systems ist. 
Im Hauptachsensystem hat ie die Komponenten — e; und — e 
Mittels der transponierten Matrix Qf, (14) erhalt man 


haben. Nun ist aber 


n* 


x (/—& 
y|=%|—e, |, 
0 0 
also 
x = — (eg cost s —e, sinT 8) , 
y =— (egsints +e, costs), 


und damit lautet die Randbedingung 


l l 
eae tds= (7) ds = 
x,(I) also ft ;] 


4. Der mathematische Inhalt des Problems. Der mathematische Inhalt des Problems besteht 
darin, das im allgemeinen inhomogene System von zwei linearen Differentialgleichungen mit nicht 
konstanten Koeffizienten, das durch (20) dargestellt wird, unter Beachtung der Randbedingungen 
(22) bis (24) aufzulésen. Da fiir jedes Stabende, also fiir s = 0 und s = / Randbedingungen vor- 
geschrieben sind, wird eine Lésung von (20) unter gleichzeitiger Erfiillung aller Randbedingungen 
nur existieren, wenn von den Parametern der Differentialgleichungen gewisse Werte angenommen 
werden; das heiBt, es liegt hier ein Eigenwertproblem vor. 

In der Regel werden wohl die Parameter a, b, c und t vorgeschrieben sein. Der fiir die Lésung 
erforderliche Eigenwert wird daher als noch verfiigbare GréBen J und q enthalten. Man kann also 
aus dem errechneten Eigenwert bei vorgeschriebener Stablange auf die Knicklast q oder bei vor- 
geschriebener Belastung q auf die Knicklange | des Stabes schlieBen. 


J (eg cost s—e, sints)ds =0, 


n 


I 
0 
(24) 
ERS ae a ao. 


5. Die Lésung der homogenen Differentialgleichungen. Erste Voraussetzung fiir die Lésung des 
eigentlichen Eigenwertproblems ist zunachst die Kenntnis der allgemeinen Lésung des zu (20) 
gehérenden homogenen Systems. Dieses lautet 


eee ee | 


$3 : (25) 
be,—(a+ b) te, + [q—s)—ar]e=0. 


Durch Einfiihrung einiger Abkiirzungen und Variablentransformationen kann man die Glei- 
chungen (25) so umformen, da man den allgemeinen Typus, zu dem sie gehéren, besser erkennt. 
Es sind dies die Abkiirzungen 


a 


3 
peetle yi ey ant (26) 


18* 
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und die Transformationen 
x= tM mae qT, 
e(8) > 4? e,(x) , (27) 
e,(s) > p(x) « 
Damit erhalt man statt (25) 
@, — (¢ + 1) el? TA-3 ce, + [x 4+ T?21-48 & *— jee = 0 | (28) 
é, + (4 +1) x? TA e,+22xe, =0, 


wobei jetzt und im folgenden Punkte Ableitungen nach x bedeuten sollen. 
Diese Differentialgleichungen gehéren zu einer Klasse, deren allgemeine Form lautet 


€ +ae, + Q(x) e =90, €; —X & + y(x)e, = 90. (29) 


Dabei ist « eine Konstante; p(x), p(x), e,(x) und e,(x) sind Funktionen der unabhangigen Variablen Xe 
Dem System (29) ist eine Differentialgleichung vierter Ordnung fiir e,(x) gleichwertig: 


yP— PP +W+ yor + po) e+ Zyp—$e—VO) & + (WP—PP TP) %& =9. (30) 
Hat man eine Lésung e,(x) von (30), so erhalt man die entsprechende Funktion e,() aus 


=e P+ (p 0%) be +o eal: () 


Das System (29) oder die damit gleichwertige Differentialgleichung (30) haben eine gewisse Be- 
deutung bei technischen Problemen, die hier in bezug auf das Allgemeine erlautert werden soll: 


Fall I: Es seien y und py Konstante. Auf diesen Fall lassen sich unter anderen die Differential- 
gleichungen folgender Probleme zuriickfiihren: 


I, 1) Der Spannungsfrei verwundene, von einer Einzelkraft an den Stabenden belastete Druck- 
stabt. 

I, 2) Die tordierte und von einer Einzelkraft an den Stabenden belastete Welle’. 

I, 3) Das Foucaultsche Pendel*. 

I, 4) Kritische Drehzahlen fiir die einfach besetzte Welle*. 


In diesem Fall I baut sich das Fundamentalsystem der Lésungen aus trigonometrischen Funk- 
tionen auf, wie allgemein bekannt ist. 


Fall Il: Es seien g und wy lineare Funktionen von x. Hierzu gehéren die oben genannten Pro- 
bleme I, 1; I, 2 und I, 4; wofern man annimmt, da zusatzlich auf Stab oder Welle ein axialer 
Druck wirkt, der von einer gleichmaBig iiber die Achse verteilten Belastung (z. B. Eigengewicht) 
herriihrt. Es gehért also insbesondere auch das in dieser Arbeit zu erérternde Problem dazu. Das 
oben genannte Problem I, 3 verwandelt sich hier in das Problem von der Zwangsbewegung eines 
Punktes bei kleiner Auslenkung aus der Ruhelage auf der Flache eines Rotations-Ellipsoides unter 
der Wirkung eines Kraftfeldes, dessen senkrecht zur (e,, e,)-Ebene stehender Kraftvektor linear 
von der Variablen x abhangig ist. Die Bewegung sei auBerdem in einem gleichférmig rotierenden 
Achsenkreuz beschrieben. Fiir diesen Fall II muB8 eine allgemeine Lésung des Systems (29) noch 
aufgesucht werden. 


Fall IIT: Es seien y und y beliebige Funktionen von x. Dann ist die bei Fall II eingefiihrte 
Druckbelastung zwar kontinuierlich, aber nicht mehr gleichmafig, sondern beliebig tiber die Achse 
von Stab und Welle verteilt, und im letztgenannten Problem ist der Kraftvektor beliebig von x 
abhangig. Dieser Fall III soll nicht weiter in Betracht gezogen werden. 


Fall IV: Es kann noch der Sonderfall eintreten, daB y =y sind. Das trifft zu, wenn die Biege- 
steifigkeiten fiir beide Haupttrigheitsachsen in allen Querschnitten von Stab oder Welle einander 


1S. FuBnote 3, Seite 262. 

2 S. SuBnote 4, Seite 262. 

3 A. Sommerfeld, Vorlesungen iiber theor. Physik Bd. I, S. 162—168. Leipzig 1944, 

4 C. B. Biezeno u. R. Grammel, Technische Dynamik Bd. II, S. 177. Berlin/Gottingen/Heidelberg 1953. 
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gleich werden, bzw. wenn die Zwangsbewegung nicht auf der Flache eines Rotations-Ellipsoides, 
sondern auf einer Kugelflache erfolgt. Fiir den verwundenen Knickstab wird das Problem dann 
trivial; trotzdem soll auch dieser Sonderfall wenigstens hinsichtlich seiner mathematischen Behand- 
lung eine kurze Betrachtung erfahren. 


Fir diese Arbeit hat nur der Fall II Bedeutung, da er auf den verwundenen Knickstab mit 
gleichformig verteilter Druckbelastung zutrifft. Denn nach (28) und (29) ist 


p(x) = [x + T2148 (44 — 1)] 23, 
Y(x) =? x, (32) 
wo = — (4 + 1) x? T -13 | 


und das sind lineare Funktionen von x, bzw. es ist eine Konstante. 
Unter Verwendung von (32) geht (30) iiber in 
a ef) — ee) + (A x + Bx?) é. + (Cx — A) ey + (Ex? + Dx?—C)e, =0 (33) 


mit den Abkirzungen 


A= wl 3 2) TAA 28, 
B=+#(1+ x), 


Co, (34) 
D =+A(1 —x) T2728, 
50%, 


Ehe an eine Behandlung von (33) herangegangen wird, soll zuerst der Sonderfall a = b (gleiche 
Haupttragheitsmomente in allen Querschnitten) behandelt werden. Die hieraus hervorgehenden 
Lésungen gelten nur fiir die tordierte und gleichzeitig gedriickte Welle, fiir kritische Drehzahlen 
und nicht fur den hier zu behandelnden Knickstab; aber man wird durch die mathematische Be- 
handlung dieses Sonderfalles eine Verwandtschaft zwischen den Funktionen des Fundamental- 
systems von (33) zu schon bekannten Funktionen erkennen, was von Interesse sein dirfte. Die 
ebenfalls vorhandene Verwandtschaft mit trigonometrischen Funktionen ist schon erwahnt worden. 


Aus (26) und (32) folgt, daB fiir a = 6b gilt 
el CS UE C= = x. 


Dann ist das homogene Differentialgleichungssystem entartet zu 


égt+ae +xe,=0, é: Xe, + xe, =0, (35) 


und die Abkiirzungen (34) werden zu 
Ae Be eo 1 DO Se, 


Die gleichwertige Differentialgleichung vierter Ordnung lautet statt (33) 
x) — eB (a2 x + 2 x”) e, + (x — a?) ey + (x? — lhe, = 0. (36) 


Fiir (36) laBt sich eine geschlossene Lésung mit bekannten Funktionen angeben. Man setze an 


ie r(x) sin > x bzw. €, = 1(x) cos > ae (37) 
Geht man mit einem dieser Ans&tze (37) in die Gleichung (36) ein, so erhalt man in beiden Fallen 
als mégliche Bedingung fiir die Erfiillung von (36) die beiden getrennten Differentialgleichungen 
fir r allein 


COCs ae CLG a 3 
x ri) — 7) +Q#—9)F- («—F)i (#— ge ygpute—ajr=0. 


3 : a 1 3 5 
a — Tree t Tat +(ge— jee |r =o. 


Sie miissen beide gleichzeitig erfiillt sein. 


“ 
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Nun lassen sich aber die Gleichungen (38) umschreiben in 
7 2 3 
r+ Gs | J] r aie [~ ai Pak *| 


% a? 
if e Fa)? 
a ee: \r|—= |? +(F +47] a)! 5 


Hieraus ersieht man als weitere Méglichkeit fiir die Erfiillung von (36), daB r als Lésung der Diffe- 
rentialgleichung 


if+(+9) 


(39) 


t+ (q+ 2)r=0 (40) 
zu bestimmen ist. Diese Differentialgleichung ist bekannt. Sie laBt sich mit 
oak (41) 
auf die Gestalt 
d’r 
ae = 4 
dz? ae Aire (42) 


bringen. Hierfiir sind die Lésungen 
9 Meee ae 1/2 2 /o2  \3/2 
ry = 22.941) G Sus *| 91/3 alae 2 | 


2 a 1/2 2 [OP 3/2 
=F 910 (3 ala =) 83[5 (4+ 4 | 


: 2 (a2 3/2 
Dabei sind 1/3 bzw. S_ 1/3 Besselsche Funktionen des Argumentes Aon. + x| . Verwendet man (43) 


(43) 


fiir den Ansatz (37), so sieht man, daB der Lésungsweg nicht nur méglich und hinreichend, sondern 
auch notwendig war; denn man hat durch (37) und (43) vier linear unabhangige Partikularlésungen 
von (36) gefunden und damit die Gleichung (36) vollstandig gelést. 


Damit ist die Verwandtschaft der Funktionen des Fundamentalsystems von (33) mit den 
Besselschen Funktionen von der Ordnung = und — 2 gezeigt. 


Hat man aber e,(«) errechnet, so erhalt man e,(x) gemaB (31) aus 


1 


e,(x) = ie [ef + (a? + x) e, + eg], (44) 
was auf die beiden Lésungsgruppen 
= (| sin > Re baw. es = (a cos + Ti, 
oder . is 4 (45) 
es == in| cos > x baw. n= fn sin = 
fiihrt. 


Jetat sei die Differentialgleichung (33) behandelt. Es ist eine lineare Differentialgleichung vierter 
Ordnung mit variablen Koeffizienten. Sie hat eine auBerwesentliche Singularitat (Stelle der Be- 
stimmtheit) bei « = 0 und eine wesentliche Singularitat bei x = oo. Hinsichtlich der Singularitaten 
verhalt sich die Differentialgleichung demnach wie eine Besselsche Differentialgleichung, mit der 
sie ja auch verwandt ist. Da keine weiteren Singularitaten auftreten, lassen sich die Lésungen 
durch Potenzreihen in x darstellen, die in der ganzen Ebene um x = 0 als Mittelpunkt konvergent 
sind. 


Geht man nach der Methode von Frobenius! vor, so erhalt man vier partikulare, linear unab- 
hangige Lésungen von der Form 


t( 2) 2 I ae (46) 


1 Hort-Thoma, Die Differentialgleichungen der Technik und Physik, S. 180—186. Leipzig 1950. 
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Die méglichen Werte von t sind dabei 4, 2,1 und 0 als Wurzeln der determinierenden Gleichung 
i(t — 1) #— 2) (#— 4) = 0. (47) 


Obwohl diese Wurzeln ganzzahlige Differenzen haben, tritt hier der besondere Fall auf, daB die 
partikularen Lésungen von (33) allesamt keine Logarithmen enthalten. Das stimmt mit der schon 


bekannten Tatsache iiberein, daB die Lésungen (45) der Differentialgleichung. (36) auch alle log- 
arithmenfrei sind. 


Die sechs ersten Koeffizienten a, von (46) sind 


Go — ul"; 
Ue 
e. A 
GEE +2)” 
ae C+ Bt 48 
: a3 C+IC+26¢43 v4 


A2 
A= ED) (922) @ 3) 4)” 
(+ 2)D—2AC—2%Xt+ DAB 
CFIO+2¢+3)¢+H0+435)' J 


Sie ergeben sich aus der Rekursionsformel 


a, — 


eg — inn 
rr na = Us 
- A Cae BG 3) 
eee k= ae = 1) RS I) GS Eh) 


PNG Gt ke =3) RG E33) 
ROS ee Ey bol): “(hk Ge Ob) Oe 


Setzt man nacheinander in (48) und (46) die vier méglichen Werte fiir t ein, so erhalt man mit 
den vier Integrationskonstanten K,, K,, K,, K; als allgemeine Lésung von (33) 


ae) Ko SP sg Ot sea See S| 
+R Psy oS tA oe eee 
+ Kye {lag See a 
+Kfi—A4e—Serfa— ee eal, (49) 


Zur Abkiirzung sei der Begriff der ,,abgewandelten oder transmutierten Bessel-Funktionen“ ein- 


gefiihrt: 


Bays) = Ah — 05 | 

Bu(x) = 22 ( Soo ce r ce 
Bile) = [l— yoy — r 

B(x) = leeiheeg 


Dann kann man als allgemeine Lésung des homogenen Systems (28) bzw. der damit gleichwertigen 
Differentialgleichung vierter Ordnung (33) schreiben 


e(x) = Ss K, B(x) , (51) 
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und gemaB (31) unter Verwendung von (32) 


1/3 : 

AE) TED fe) 4 [pdx + T2479 x(1 + 3x)] by + eg}. (52) 
Zu beachten ist noch, da8 man beim Ausschreiben von e,(x) nach (51) oder in (52) fiir die Anwendung 
in der praktischen Rechnung die Konstanten A bis E gemaf (34) auszuschreiben hat. 


6. Die Eigenwerte des Knickstabes. Bei den Knickfallen 2 und 3 von Abb. | sind die Stabenden 
querkraftfrei, und die Konstanten A*, B* sind nach (21) gleich Null. Die Differentialgleichungen (20) 
werden also homogen und lassen sich ohne weiteres in die Form (28) iiberfithren. Die Lésungen des 
Problems sind daher unmittelbar durch (51) und (52) gegeben. 

Bei den Knickfallen 1 und 4 sind die Stabenden dagegen nicht querkraftfrei; 4*, B* sind von 
Null verschieden. Die Lésungen (51) und (52) des homogenen Problems kénnen jetzt nicht unmittel- 
bar verwendet werden; denn die Differentialgleichungen (20) sind inhomogen und gehen mit den 
Abkiirzungen und Transformationen (26) und (27) statt in die Form (28) tiber in 
é, —(x% + 1) xl? TA-'3 ec, + [x + T24—-28 (4-1 — 1)] He, 

= — A* sin [T (1 —x% 2-3 x — J T?] + B* cos [T (2 — 24d — A 1], 
ét+(@4+)x#? TA 8 e+ xe, 

= 4-2 A* sin [T (1 —x A-V3 x — A T?)] + x1? B* cos [T (1 —xA—3 x —14T9)]. 


(53) 


Die allgemeinen Lésungen e, und e, dieser Gleichungen (53) kann man aber vermittels der 
bekannten homogenen Lésungen (51) und (52) allein durch Quadraturen (z. B. nach der Methode 
der Variation der Konstanten) finden. Da dies keine prinzipiellen Schwierigkeiten bietet, soll es 
hier nicht weiter verfolgt werden und es sollen auch die allgemeinen Liésungen des inhomogenen 
Falles vermége (51), (52) grundsatzlich als bekannt angesehen werden. 


a) Der homogene Fall. Gegeben sind 
x= = = Verhaltnis der Haupttragheitsmomente des Stabes, 


T =11 = Gesamtverwindung des Stabes; 
gesucht ist 


3 
A= ae als Eigenwert des Problems; 


vorgeschrieben sind gema8 Ziff. 3 durchweg homogene Randbedingungen: 
R, |e,» eg]. = 90; 
R,[eg, e,] = 0 


Ri len» ee] ae fai are 
Riles e)= 6 


| fae he 
(54) 


Dabei sind die Randbedingungen R;[e,, é3] (i = 1,2; 8, y =é,7) lineare Verkniipfungen von e, 
und deren Ableitungen. e 

Vermittels der Beziehungen (27) kann man die Randbedingungen (54) leicht umschreiben in 
solche fiir e,, e, und x: 


Ry[eo, e)=0, 
Role), é] = 0 
Rifles, é)=0, 
R[e,, é.] = 0 


fir x =x—1 413 (1—{-1 T2), 
(55) 
fiir x= — 9-1-4372, 


Wendet man die Randbedingungen in der Form (55) auf die Lésungen (51) und (52) an, so 
erhalt man im homogenen Fall bei vier Integrationskonstanten K, bis Kg ein vierzeiliges lineares 
homogenes Gleichungssystem 


L(K) = a9 Ky + a9; Ky + A. K, + ag, Kz, 
L,(K) = a4) Ky + a,, K, + a,. K, + a,, Kj, 
L,(K) = dg) Ky + 5, Ky + a9 Ky + agg K;, 
L,(K) = ago Ky + a3, Ky + 39 Ky + a35 Ky 


(56) 
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zur Bestimmung der K; (i = 0,...,3). Die Koeffizienten G3 (U = 0;...5.3; 6=0,..., 3) sind 


dabei Ausdriicke, die sich aus den GréBen x, 2 und T zusammensetzen. 


In bekannter Weise muf fiir eine eindeutige Bestimmung der K, die Matrix des Systems (56) 
singular, also die Determinante gleich Null sein: 


Mo Ay AQ Ags | 

| 

1449 %1 Ge Gz) 
det b= | 0 (57) 


| 220 Go, Age G3 


| G39 %3 1 Age Ags 
Diese Determinante ist aber ein Ausdruck, der, wie seine Elemente, die GréBen x, A und T enthalt: 
detikc== Ox A. T).= 0. (58) 


Hierbei sind x und T gegebene GréBen, so da (58) eine Bestimmungsgleichung fiir den Eigenwert A 
darstellt, aus welcher man diese gesuchte GréBe ausrechnen kann. 


b) Der inhomogene Fall. Er unterscheidet sich von dem vorhergehenden Fall a) dadurch, 
daB jetzt gemaB (24) noch zwei Randbedingungen 


R,le, et] = 0, R,[e,, &,T] = 0 


fiir s = 1 bzw. x = — x—14—?/3 T? hinzukommen. Ferner sind die Randbedingungen nicht auf die 
Lésungen (51) und (52), sondern auf die mit ihnen ermittelten inhomogenen Lisungen anzuwenden, 
Dann hat man, wie es sein muB, ein sechszeiliges lineares, homogenes Gleichungssystem zur Be. 


stimmung der sechs Integrationskonstanten K, bis K, und A*, Be: 


Alles tbrige lauft, wenn auch mit mehr Arbeitsaufwand, wie vorher: aus der gleich Null ge- 
setzten Determinante des linearen Gleichungssystems folgt eine Bestimmungsgleichung fir den 
Kigenwert A. 


7. Die numerische Rechnung. Fiir die zahlenmafige Ermittlung von Eigenwerten werden die 
gem4B (51) und (52) in Reihenentwicklung vorliegenden Lésungen e, und e, wohl unhandlich sein. 
Der in Ziff. 6 geschilderte Lésungsweg fir die Eigenwerte 4 diirfte daher wohl nur grundsatzliche 
Bedeutung haben. 


Fur die wirkliche Durchfiithrung der Rechnung soll daher hier ein anderer Weg beschrieben 
werden. Dabei soll der Einfachheit halber der homogene Knickfall 2 von Abb. 1 zu Grunde gelegt 
werden. Es wird nicht schwer sein, in einem anderen und auch inhomogenen 
Knickfall entsprechend vorzugehen. 


Fiir den Knickfall 2 (Abb. 5) gelten das homogene Gleichungssystem (25) 
und die Randbedingungen (23) fiir beide Stabenden. 


Die Gleichungen (25) werden durch 


Gee eT od 
‘ (59) 
Tied == 
a l 
iibergefiihrt in 
é, + (1+) Teg + AL —o)—xT]e,=0, 
dees a Aes are (60) 
es — Le A —o) — 1" « Fer —0, 


4 
wo jetzt und im folgenden Punkte Ableitungen nach o bedeuten. Die Randbedingungen (23) gehen 


dadurch tiber in 
(0) ae fiir = 0, (61) 
e(o) = 0 lie 
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Auf die Differentialgleichungen in der Form (60) und die Randbedingungen in der Form (61) wird 
die Differenzenmethode angewendet. Denkt man sich den Stab in n Intervalle unterteilt, so erhalt 
man 


aa ne = 0, 
” ” n Ton 


gt Diet Ae ae if (lin) (ee — era) 21 on lat Hh 


(62) 


ln *e,5= OF 
& Din 4 al 7] é 


0, 


Feta ac pf) 


was ein finites, lineares Gleichungssystem fiir die Ordinaten der Naherungen von ég, @, darstellt; 
i lauft von 0 bis n. 


Durch Ausrechnen der beiden oberen Gleichungen (62) erhalt man die Rekursionsformeln 


oa Dee ape, Uae ag SV Si Aas \ear 4uW 4nd. ee 

Uh at gzaaaargia Lee +( nZ — az) #3 Z nZ *) &n? (a) 
ans qgee ons 1 24 

eit] — — gi—1_“_* gi—1__/—______ {| ek —(—_ —— 51] ¢,. 

§ ZOE he SD Z nZ 5 nZ n? Z a 


Dabei sind die folgenden Abkiirzungen verwendet worden: 


a ae 

U = 7T*(1 +2?—4n?x, 

V=A— T?—2n'*x, (64) 
W =j)—x T?—27n’, 

Z=T1+xy?+4n?x. 


Das Eigenwert- und damit Randwertproblem wird nunmehr auf ein Anfangswertproblem zu- 
rickgefiihrt. Das Differentialgleichungssystem (60) enthalt den Eigenwert / als Parameter, folglich 
werden die Lésungen von (60), bekannterweise, Funktionen dieses Parameters sein. Durch Variieren 
des Parameters A wird es dann méglich sein, gerade die Lésungen von (60) zu finden, die die Rand- 
bedingungen (61) erfiillen; das sind die Eigenfunktionen des Problems; der zugehérige Wert A ist 
gerade der gesuchte Eigenwert. 


In der Durchfiihrung sieht das so aus: Man nimmt einen bestimmten Zahlenwert fiir 4 an; 
x und T sind gegebene Zahlenwerte der Aufgabe. Damit errechnet man die GréBen (64). Jetzt 
verwendet man als Anfangswerte die Randbedingungen e? = e)? = 0. Da die Gleichungen (62) 
homogen sind, kann noch der Freiwert e 1 = 1 angenommen werden. e, + bleibt als unbestimmte 
GréBe in der Rechnung. 


Nun kann man die Rekursionsformeln (63) anwenden und schrittweise alle Ordinaten e%, e} 
von t = 0 — n aufbauen. 


Wendet man dann z. B. auf e’ die vorgeschriebene Randbedingung 
Cea 0; also ee e—a() (65) 


an, so folgt daraus eine Bestimmungsgleichung fiir die bisher unbestimmte GréBe e7 1. Hat man 
ihren Wert errechnet, so kann man ihn in den Ausdruck fiir ez einsetzen. Wird dieser Wert zu Null, 
wie es dic letzte Randbedingung am Stabende vorschreibt, so war der Wert von A gerade der ge- 


winschte Kigenwert. Wird er dagegen nicht zu Null, so tragt man ihn als Fehler f in einem f, /- 
Achsenkreuz auf. 


Fihrt man dies fiir eine Reihe von angenommenen /-Werten durch, so erhalt man eine Fehler- 


funktion f(A). Die Nullstellen dieser Funktion geben die Eigenwerte 1 mit einer dem Differenzen- 
verfahren zukommenden Genauigkeit. 


Fir die praktische Durchfithrung der Rechnung ist es vorteilhaft, die Anwendung der Rekur- 
sionsformel (63) als wiederholte Matrizenmultiplikation aufzufassen!. 


1 R. Zurmiihl, Matrizen, S. 382—384, Berlin/Géttingen/Heidelberg 1958. 
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Es werden eingefithrt der Zustandsvektor: 
e—1 
we 
— oa (66) 
ee 
et 
und die Punktmatrix 
0 0 1 0 
0 0 0 i 
= 67 
’ P31 Ps2e Pss Psa 8 
Pai Pa2 Pas Paa 
mit den Elementen 
U 2nznS 
P31 = 7-2 Pal Se ga ] 
2nS U 
Lae a> gre Pita sy 3 
4V Ah. SV SAe (68) 
Eee 3? Pa GZ ae 
SW Sin 4nW Anda. 
Pee 7 Pie oe 


wobei die Elemente p33, P34, Pag Pag Sich mit der von Schritt zu Schritt veranderlichen GréBe i 


andern. 


Wendet man die Vektoren 3' bzw. Matrizen ‘3; auf die Rekursionsformel (63) an, so erhalt man 


statt dieser 


aitt 


ier 


(69) 


_ Lauft i von 0 bis n — 1, so erhalt man als Resultat der wiederholten Matrizenmultiplikation 


i Br Bn—s ome $B %, Bo oo 


(70) 


Man kann das Produkt %, 1 B,2... Bo $8, $y der Punktmatrizen noch zusammenfassen zur 


Ubertragungsmatrix 

Ura Mie slag. ta 

ul oat Ups nen We 

Us) goes; seale 

Wane) nip = Urigy, Ur 

Der Anfangsvektor 

e;- 
p=(% 
ce 
en 


== gba eo one Be ¥, Bo - 


(71) 


_ kann beziiglich seiner Elemente naher bestimmt werden. Auf Grund der Randbedingungen ist 


wegen der Homogenitat der Gleichungen (62) kann 


ae Oe 
ee =e = 90; 


_stimmt. Damit kann man fiir 3° schreiben 


et = 
1 1 
rg eg ty oat 
0 0 ae 
0 0 


1 gesetzt werden. Nur e, 1! bleibt unbe- 


(72) 
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Auch bei 


ee 
ee (73) 
0 
0 
Verwendet man (71), (72) und (73), so wird aus (70) 
er i) eee ee 1 
en a | Vero Ven Mas Me en ; (74) 
0 Us, 2 ep es Me3 uke 0 
0 gi Ua Cage es 0 


Beachtet man bei (74) die durch die Randbedingungen am Stabende vorgeschriebenen beiden 
letzten Zeilen, so erhalt man 


Us; + Ugg e, = 0, Ug, + Ug2e, = 9. (75) 
Hieraus folgt a 
ert = — Mal ; (76) 
” Te 


Sollte man den zur Berechnung der p,, nach (68) und damit der u,, nach (71) verwendeten 
Betrag des Parameters J gerade so passend gewahlt haben, daB auch noch 


Us, + Usee," = 0 (77) 


ist, so hat man mit diesem Wert von / einen Naherungswert fiir den Kigenwert gefunden und kann 
die bis auf einen gemeinsamen Faktor bestimmten Ordinaten der Eigenfunktionen auf Grund 
von (69) ausrechnen, weil man ja nun die Rechnung mit einem Vektor 3° beginnen kann, bei dem 
in (72) fiir e7 1 der gemaB (76) bekannte Zahlenwert eingesetzt werden kann. 


Diese allgemeinen Ausfiihrungen seien durch ein Zahlenbeispiel erlautert. Es sei ein Stab vor- 
gelegt mit x = b/a = 5. Fir die Unterteilung des Stabes sei n = 10 angenommen. Die Rechnung 
sei fiir verschiedene Gesamtverwindungen T des Stabes durchgefiihrt, und fiir jedes bestimmte T 
erhalt man einen anderen EKigenwert 47. Fiir T = 0 (das ist der verwundene, gewohnliche Knick- 
stab unter gleichférmig verteilter Belastung, fiir den die Eigenwerte mittels Besselscher Funktionen 
errechnet werden! hat man als erste beide Eigenwerte 


f= 18.05 eI 818: 


Mit allen Eigenwerten Af, wird das Verhaltnis 


=e | (78) 


gebildet und dieses dann als Funktion von T aufgetragen. Dieses Verhaltnis gibt an, um wieviel der 
Eigenwert gréBer als der des unverwundenen Stabes ist; da stets k > 1 ist, sieht man, daB® durch 
das Verwinden des Stabes die Tragfahigkeit oft erheblich heraufgesetzt wird. 


Fir T sollen die Werte 0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°, 360° angenommen werden. 
Ist z. B. T = 180° so hat man fiir die Ausdriicke (64) und daraus fiir die Elemente Pa,p der 
Punktmatrizen ‘f; (68) in Abhangigkeit von 4 die in der Tabelle 1 angegebenen Werte. 


_ 1 Fr. A, Willers, Z. angew. Math. Mech, 21 (1941) S. 43. 
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Tabelle 1. T=180°; x=5; u = 10 
S = 37,699; U = — 1644,680; Z = 2355,320 
Pai P32 Ps3 Ps Pai Ps Pas Pas 
meat oe 
| 

i * 7 U peta A) 4h SWE ines | U | SK Sz |. 4uW 4x) 
5 ef sot if es coeds sid = eles 
Zz he Zz nZ nZ nz Zz Z n Z nL, Z, nZ 
20 —989,870 —229,350 | —0,698 | —0,320| 1,681+0,003 i 0,367-+-0,003 i 1,600 | —0,698 | —1,584—0,0037| 1,948+-0,017i 
25 —984,870 —224,350 | —0,698 | —0,320| 1,673+0,004 i 0,359 + 0,004 7 1,600 | —0,698 | —1,576—0,004i/} 1,905+0,0217 
30 —979,870 —219.350 | —0,698 | —0,320| 1,664+-0,005 i 0,3514-0,005 7 1,600 | —0,698 | —1,568—0,005i| 1,863-+0,025 i 
35 —974,870 —214,350 | —0,698 | —0,320| 1,656+0,006i 0,343+-0,006 i 1,600 | —0,698 | —1,560—0,006i| 1,820+0,030i 
40 —969,870 —209,350 | —0,698 | —0,320| 1,647+.0,007i 0,335 +-0,006 i 1,600 | —0,698 | —1,552—0,006i/} 1,778+0,034 i 
45 —964,870 —204,350 | —0,698 | —0,320}| 1,639-+0,008 i 0.327+-0,007 i 1,600 | —0,698 | —1,544—0,0077| 1,735+-0,038 i 
50 —959,870 —199,350 | —0,698 | —0,320| 1,630-+0,008 i 0,319+.0,008 i 1,600 | —0,698 | —1,536—0,008i| 1,693+-0,042 i 
55 —954,870 —194,350 | —0,698 | —0,320! 1,622+0,009 i 0,3114-0,009 i 1,600 | —0,698 | —1,528—0,009i| 1,650+0,047 i 
60 —949,870 —189,350 | —0,698 | —0,320| 1,613+0,010i 0,303 +-0,010 i 1,600 | —0,698 | —1,520—0,010i|} 1,6084+0,051i 
65 —944,870 —184,350 | —0,698 | —0,320) 1,605+-0,011i 0,2954+-0,010 i 1,600 | —0,698 | —1,512—0,010i| 1,565+0,055 i 
70 —939,870 —179,350 | —0,698 | —0,320| 1,596+0,012i 0,287-+0,011 i 1,600 -| —0,698 | —1,504—0,011i| 1,523+-0,059 i 


Tabelle 2. A =50; T = 180° 


—1 -1 
1 en it en 
1 0 0 0 1 0 2,702 ~ —0,660 | 
0 1 3 oy 0 0 0 1 5,432 1,244 | 3° = 
1.0 0 ut ° | —0,698 —0,320 1,670 0,359] 3,457 —0,174| Bs: 4° 
He ekO) 0 1,600 —0,698 —1,576 1,903} 4,266 1,349 
0 0 4 0 0 0 ; lr 20 0 1 o | 3,457 —0,174 | ; 
0 0 1 0 0 = vie. 0 0 1 4,266 1,349 | 3° = 
Bo —0,698 —0,320 1,630 0,319 | —0,698 —0,320 | Wo: 3° Be | —0,698 —0,320 1,678 0,367 | 3,742 0,266 | e- 3° 
1,600 —0,698 —1,536 1,693 1,600 —0,698 | 1,600 —0,698 —1,584 1,945 | 3.353. 0,975 
0 0 1 0 | —0,698 —0,320 0 0 1 0 3,742 0,266 
oy ) 0 0 1 1,600 —0,698 | 4° = B 0 0 0 1 3,353. 0,975 | 4° = 
* | —0,698 —0,320 1,638 0,327 | —0,620 —0,752| 1: 3% 7 | 0,698 —0,320 1,686 0,375 | 3,788 0,504 | ‘Br- 37 
1,600 —0,698 —1,544 1,735 3,854 —0,717 1,600 —0,698 —1,592 1,987 | 3,259 0,294 
on 0 1 0 —0,620 —0,752 isaac 0 1 0 3,788 0,504 | 
oy 0 0 0 1 3,854 —0,717 | 3° = | 20 0 0 To) 3,259 03294 ig 
* | —0,698 —0,320 1,646 0,335 0,246 —1,031 | B2- 3 * | —0,698 —0,320 1,694 0,383 | 3,980 0,469 | Bs: 3° 
| 
1,600 —0,698 —1,552 1,777 Earn Osi p) | 1,600 —0,698 —1,600 2,029 4,199 —0,465 | 
: 0 0 1 0 0,246 —1,031 0 0 1 0 3,980 0,469 
oN 0 0 0 1 5,577 —0,132 | 34 = 0 0 0 1 4,199 —0,465 
* | —0,698 —0,320 1,654 0,343 1,519 —0,996 | Bs 3° > Usy Usage lego 
- | 1,600 —0,698 —1,560 1,819 6,079 0,666 —0,698 —0,320 1,702 0,391 | 4,729 0,171 | By 3° 
; Ug Ugo 
) 0 1 0 1,519 —0,996 1,600 —0,698 —1,608 2,071 | 6,082 —1,116 
mS 0 0 0 1 6,079 0,666 | 3° = 
* | —0,698 —0,320 1,662 0,351 2,702 —0,660 | Ya: 3* ey Un 6,082 
1,600 —0,698 —1,568 1,861 5,432 1,244 &, a + 7376 = 5450 
42 ’ 


f = Ug, + Use et = 4,729 + 0,171 - 5,450 = 5,661 


Der erste Zustandsvektor, mit dem die Rechnung begonnen werden kann, ist gemaB (72) vor- 
_ geschrieben. 
Samtliche Elemente p,, der Punktmatrizen 8; kénnen aus der Tabelle 1 herausgelesen werden. 
Die Rechnung nach der Vorschrift (69) wird in Tabellenform durchgefiihrt. Dies sei einmal fiir 
den Fall T = 180°; 2 = 50 in der Tabelle 2 vorgefiihrt. 
Aus der Tabelle 2 liest man : 


- ab, also nach (76) 
265082 


e, — 1,116 = 5,450 ° 


; Statt (77) erhalt man nicht Null, sondern den Fehler 


f(A) =f (50) = ug — Ug, 


9 
U4g2 


der mit den Werten in der Tabelle 2 den Zahlenwert 
(50) = 4,729 + 5,450 - 0,171 = 5,661 


. 
— 
= 
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hat. Fiihrt man nun ganz entsprechende Rechnungen fiir andere Werte von A durch, so erhalt man 


die Fehlerkurve Abb. 6. 
Aus der Fehlerkurve liest man als Abszissen der ersten beiden Nullstellen die ersten beiden i| 
Eigenwerte ab zu 


7180 30,0; 218 ww 64,8. 


Nach (78) gibt das 
30,0 


64,8 
ean ee ie 
18,9 3788 


18,9 


piso — 158, ks = 


Fiihrt man die entsprechenden Rechnungen fir x = 5, n = 10 und die tbrigen Werte von | 


0° < T < 360° durch und tragt die zugehérigen k7-Werte als Funktion von T auf, so erhalt man |f 


Abb. 7. 
(A) 
40 
7 
30 5 
20\- 
Y 
10\- 
#(50) 
peee are Cow l A. 
2 30 35 WO 45 50 55 60 65 3 
| 
-10 bo 
2 
ef) 
oe gle {ae cd ae 
O° 45° 90° 135° 180° 205° 210° 35° 360° 
Abb. 6. Fehlerkurve f(A) fiir T = 180°. Abb. 7. 
tut 
GE) en 


WS ®EA DY 


Abb. 8, Eigenfunktionen fiir T = 180°; A, = 30. 


In gleicher Weise kénnte man fiir andere Knickfalle und andere x-Werte vorgehen. Man erhielte 
dadurch eine ganze Schar von Kurven. Ihre Verwendung in der Praxis geschahe folgendermaBen. 
Dem Knickfall und dem ~-Wert entsprechend sucht man die zustandige k?-Kurve Ee . 
Der Gesamtverwindung T entsprechend liest man aus ihr k” heraus. Hiermit bildet man 


qi 
i AG ep Loy 


rer A eee Tp eS 
nab _ Av aN { 
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Die Biegesteifigkeit a des Stabes wird ebenfalls eine gegebene GréBe sein, so daB man bei bekannter 
-Stablange | fiir die kritische Last q findet 


_ ak; - 18,9 


' (78) 


was das Ziel der ganzen Rechnung ist. 


Tabelle 3. A = 30; T = 180° 


eae | ls 
| 1 5 n : re | fi t en | 
1 0 0 0 1 Olen 3:719) 1,034 
0 Lea Rees peal ii.0 0 0 1/1 5,268 1,307) 3°= 
0 0 | & ° | —0,699 —0,320 1,690 0,370| 4,597 —0,531| Bs*3 
0 0 1,600 —0,699 —1,593 1,993) 3,636 1,914 
aeO 0 1 0 0 peal 0 0 1 0 4,597 —0,531 
Bo 0 0 0 1 0 Osea R,| 9° 0 0 1 3,636 1,914) 37 = 
/—0,699 —0,320 1,665 0,350 | —0,699 —0,320 | By - 3° * | —0,699 —0,320 1,695 0,374| 4,866 0,120 | Be- 3° 
1,600 —0,699 —1,568 1,863 1,600 —0,699 | 1,600 —0,699 —1,598 2,019 | 2,263 2,145 | 
ie KE 0 1 0 | —0,699 —0,320 | 0 0 1 0 | 4,866 0,120 
Rg, | 0 0 0 1 1,600. —0,699 | 32 = oy 0 0 0 1 2,263 214504 — 
—0,699 —0,320 1,670 0,354 | —0,601 —0,782) i - 3! 71! —0,699 —0,320 1,700 0,378 | 4,751 0,773 | 3B, - 3? 
1,600 —0,699 —1,573 1,889 | 4,122 —0,817 1,600 —0,699 —1,603 2,045 | 1,641 2,007 | 
[S50 0 1 0 —0,601 —0,782 | | 0 0 1 Cee veAn rol 0,773 | 
qg,} 0 0 0 1 4,122 —0,817 | 3° = sei © 0 0 Lees) 1,641; 552500 Tes 
| —0,699 —0,320 1,675 0,358 0,446 —1,155 | B.- 3? * | —0,699 —0,320 1,705 0,382 | 4,602 1,314 | Bs: 38 
1,600 —0,699 —1,578 1,915 6,605 —0,354 | 1,600 —0,699 —1,608 2,071 | 1,963 1,606 
0 0 1 0 0,446 —1,155 ) 0 1 0 4,602 1,314 
R 0 0 0 1 | 6,605 -—0,354 | 3* = 0 0 0 1 1,963 1,606 
—0,699 —0,320 1,680 0,362 2,241 —1,260| 33° 3° I, Us Usa 310 = 
1,600 —0,699 —1,583 1,941 | 6,472 0,461 | —0,699 —0,320 1,710 0,386 | 4,781 1,684 | %B, - 3° 
0 0 1 0 2,241 —1,260 | 1,600 —0,699 —1,613 2,097 3.148 1082 
¥ 0 0 0 1 6,472. 0,461] = 
* | —0,699 —0,320 1,685 0,366 3,719 —1,034| Bu: 3* Es 3.148 
1,600 —0,699 —1,588 1,967 5,268 1,307 | oe =e : 2,909 
1 Uys 1,082 
if Sou, 5, * oe = 4,781 — 1,684. 2,909 = — 0,118 
Tabelle 4. Eigenfunktionen fiir 2 = 30; T = 180°; e;1 = — 2,909 
che ret é 
t=— |] 1,000 + 0,000 - 2,909 = 1,000 i=— 1 0,000 — 1,000 - 2,909 = | —2,909 
0 0,000 + 0,000 - 2.909 = 0,000 0 0,000 + 0,000 - 2,909 = 0,000 
1 —0,699 + 0,320 - 2,909 = 0,238 1 1,600 + 0,699 - 2,909 = 3,633 
2 —0,601 + 0,782 - 2,909 = 1,674 4 4,122 + 0,817 - 2.909 = 6,499 
3 0,446 + 1,155 - 2,909 = 3,806 3 6,605 + 0,354 - 2,909 = 7,035 
4 2,241 + 1,260 - 2,909 = | 5,906 4 6,472 — 0,461 - 2.909 = olor 
5) 3,719 + 1,034 - 2,909 = 6,727 5 5,268 — 1,307 - 2,909 = 1,466 
6 4,597 + 0,531 - 2,909 = 6,142 6 3,636 — 1,914 - 2,909 = | —1,932 
ue 4,866 — 0,120 - 2.909 = 4,517 a 2,263 — 2,145 - 2.909 =e 3d 
8 4,751 — 0,773 + 2,909 = 2,002 8 1,641 — 2,007 - 2,909 = | —4,197 
9 4,602 — 1,314 - 2. 909 = = 0,780 9 1,963 — 1,606 - 2,909 = | —2,709 
10 4,781 — 1,684 - 2,909 = | —0,118 10 3,148 — 1,082 - 2,909 = 0,000 


Von Interesse ist noch das Aussehen der Eigenfunktionen. Nimmt man / = 30 als ersten 
Kigenwert /18° an, so kann man mit ihm die Tabelle 3 aufstellen. Wie vorher errechnet man 
= ae U41 3,148 


n Stet. 1082 ee 


é 


und mit diesem Wert kann man nach (69) die Ordinaten der beiden Kigenfunktionen e, und eg 
ausrechnen, wie dies Tabelle 4 zeigt. Das graphische Bild der Eigenfunktionen ist in Abb. 8 dar- 
- gestellt. 


(Eingegangen am 26. Oktober 1959). 
Anschrift des Verfassers: Dipl.-Math. Horst Leipholz, Stuttgart-Untertiirkheim, Fiechtnerstr. 91. 
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Die Eigenfrequenzen verwundener Schaufeln* 
Von B. Jager 


Diese Untersuchung ist Herrn Professor Dr.-Ing. E.h. K. von Sanden zum 75. Geburtstag gewidmet i 


1. Einfiihrung. Im folgenden wird ein Verfahren zur Berechnung etwa der untersten vier Eigen- 
frequenzen von einseitig (elastisch oder starr) eingespannten und von gelenkig befestigten Schaufeln 
beschrieben. Torsionsschwingungen sind in der Untersuchung nicht enthalten. Es soll die Methode 
der Ubertragungsmatrizen)?? angewendet werden. Das Problem ist fiir die Durchfihrung der 
Rechnung auf einem digitalen Rechenautomaten dargestellt. Durch die Anwendung von Uber- 
tragungsmatrizen unterscheidet sich diese Untersuchung von anderen Veréffentlichungen*®®”®. 


Der Rechnung liegt die Annahme zugrunde, da die Schaufel sich wie ein Biegestab verhalte. 
Die Querschnitte quer zur Langsachse sollen in sich starr bleiben. Die Schaufel sei also durch eine 
im nicht deformierten Zustand gerade Linie ersetzt, welche mit Masse und Elastizitat belegt ist, 
wobei aber die Biegeweichheit von der Richtung der Biegung abhangt, da die Schaufelquerschnitte 
,unrund“ sind. Das Resultat der Rechnung wird also umso genauer mit der Wirklichkeit iiberein- 
stimmen, je schlanker die Schaufel ist. ErfahrungsgemaB geniigt diese Annahme fiir Schaufeln 
von Axialverdichtern mit einem Verhaltnis der Schaufelhéhe (in radialer Richtung) zur mittleren 
Profillange von etwa 4,5: 1 fiir die untersten vier Schwingungsgerade. 


Wollte man die Eigenfrequenzen noch weiterer Schwingungsgrade berechnen, so miiBte die 
Schaufel als Platte bzw. Schale betrachtet werden, da dann die Querschnitte nicht mehr als in sich 
starr angenommen werden kénnen. Bei extrem kurzen Schaufeln mit einem Verhaltnis von etwa 
1:1 treten aus diesem Grunde schon bei der Grundschwingung Abweichungen auf, vor allem bei 
stark gewélbten diinnen Profilen. Bei Verdichterschaufeln ist die Rechnung genauer als bei Tur- 
binenschaufeln, da diese im allgemeinen starker gewélbt sind. 


Da es also keinen Sinn hat, mit der Theorie der Biegung diinner Stabe mehr als die untersten 
vier Eigenfrequenzen erfassen zu wollen, kann man mit hinreichender Genauigkeit die Schaufel 
in eine Anzahl Abschnitte unterteilen, welche aus masselosen elastischen Teilstiicken und aus Punkt- 
massen bestehen. Dies entspricht auch insofern den Gegebenheiten der Praxis, als von der aero- 
dynamischen Berechnung her die Schaufel in einzelnen Profilschnitten gegeben ist. Man wird 
also an diesen Stellen die Daten der Schaufel (Flache, gréBtes und kleinstes Tragheitsmoment, Lage 
der Haupttragheitsachsen) aus den Profilschnitten berechnen. Damit wird die Schaufel zunachst 
als ein stufenartig unterteiltes Gebilde aufgefaBt, bei dem sich die geometrischen Daten (auch die 
Verwindung) sprunghaft jeweils an der Trennstelle zwischen zwei Abschnitten andern. Im all- 
gemeinen wird eine Unterteilung in zehn Abschnitte ausreichen. 


2. Die Schaufel mit starrer Einspannung am Schaufelfuf. a) Allgemeines. Fiir die Schaufel 
sei das Koordinatensystem der Abb. | zugrunde gelegt. Darin bedeutet z die Langsrichtung (Achse 
der Schaufel in radialer Richtung), y die Richtung der Welle der Maschine und x die Umfangs- 
richtung. Die Schwerpunkte aller Schaufelschnitte senkrecht zu z sollen in nicht deformiertem 
Zustand der Schaufel auf der z-Achse liegen. Eine Abweichung der z-Achse von der radialen Rich- 
tung infolge des Ausgleichs? der aerodynamischen Krafte kann hier ohne Beachtung bleiben. Die | 


* Die Untersuchung wurde bei der Firma Daimler-Benz AG in Stuttgart durchgefiihrt. Das Programm 
a oe ee enciechner IBM 650 wurde von Herrn Dipl.-Ing. R. Lorenz von der gleichen Firma 
auigestelit. 

1 H. Fuhrke, Ing.-Arch. 23 (1955) S. 329. 

K. Marguerre, J. Math. Phys. 35 (1956) S, 28. ° 

E. Pestel, G. Schumpich, S. Spierig, VDI-Berichte 35 (1959). 

E. Maier, Ing.-Arch. 11 (1940) S. 73. 

A. Mendelson, S. Gendler, NACA-TN 2300, Washington 1951. 

D. D. Rosard, J. Appl. Mech. 20 (1953) S. 241. 

F, Njordson, Acta Polytechnica, Mechanical Engineering, Series 3, Nr. 3 (1954). 
J. C. Houbolt, G. W. Brooks, NACA-TN 3905, Washington 1957. 

* Es ist tiblich, die Schaufeln soviel gegen die radiale Richtung geneigt einzubauen, daB sie bei der Héchst- 
drehzahl frei von Biegebeanspruchungen sind, in dem die Gaskrafte und die Fliehkrafte sich die Waage halten. 


On Oa er &O WD 
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Koordinaten £ und 7 seien die Haupttragheitsachsen eines Schaufelschnittes; y sei der Winkel 


_ zwischen den Achsen x und & bzw. y und n. Die Schaufel sei in n gleiche Abschnitte von der Lange 


Az unterteilt. Die Punktmasse liege am Ende 
jedes Abschnitts. Dann wird beim untersten 
Abschnitt die halbe Masse an der Einspan- 
nung vernachlassigt (Abb. 2). Beim obersten 
Abschnitt wird nur die halbe Masse eingesetzt. 
Die Flachen S; werden an den Stellen z, ge- 
nommen. Fir die Tragheitsmomente I,; und I te 
sind geeignete Mittelwerte einzusetzen (naheres 
siehe Ziff. 6). Zur Berechnung der Schwingungen 
von verwundenen Staben! sind Matrizen von der 
achten Ordnung notwendig entsprechend den zweimal vier Vektorkomponenten (Verschiebung, Nei- 
gung, Moment M und Querkraft Q in zwei zueinander senkrechten Richtungen) eines Zustandsvektors 


e 


! 
| 
! 
{ 


VILA OT 


Abb. 1. Koordinatensystem. Abb. 2, Einteilung der Schaufel. 


a Q, u 
Mit Ricksicht auf den spater betrachteten Einflu8 der Fliehkraft, welcher im Koordinaten- 
system (x, y) definiert ist, da der Einflu8 in den beiden Richtungen verschieden gro8 ist, wird der 
obige Vektor im Koordinatensystem (x, y) angeschrieben. Nach der Methode der Ubertragungs- 
matrizen kann man von dem Biegezustand an der Stelle i — 1, beschrieben durch den Vektor Deseys 


- auf den Biegezustand an der Stelle i mit dem Vektor 1); durch Multiplikation mit der Abschnitts- 
matrix YY; ; , tibergehen gemaB der Gleichung 


ji = Miia Yi-1- (2) 

Die Abschnittsmatrix stellt den Ubergang iiber ein masseloses elastisches Stiick mit der Elastizi- 

tatsmatrix & und den Ubergang tuber eine Punktmasse mit der Massenmatrix J¢ dar. Die Ver- 

windungsmatrix % beriicksichtigt die Lage der Haupttragheitsachsen des Schaufelquerschnitts. 
Mit diesen Matrizen lat sich fur die Abschnittsmatrix schreiben 

eee ES OS Te (3) 


b) Die Elastizitatsmatrix {. Diese sei in den jeweiligen Haupttragheitsachsen (€, 7) des 


 Querschnitts, welcher fiir einen Abschnitt maBgebend ist, definiert. Bei einer Biegung, welche 


in den Haupttragheitsachsen erfaft ist, verschwinden die Koppelungselemente der Matrix, was fiir 
die Elemente des rechten oberen und des linken unteren Viertels der Matrix zutrifft. lm linken 
oberen Viertel stehen nach Fuhrke? die Elemente fiir eine Biegung in einer Richtung, namlich 
der Achse 7 mit dem Tragheitsmoment I, (um die §-Achse), im rechten unteren Viertel die Elemente 
fiir eine Biegung in der Richtung € mit dem Tragheitsmoment J,. Die Matrix lautet nun 


e Az Az?/2 E I, Az/6 E I; | 


eee As | : 
0 0 0 ae 
L= : (4) 
bicreAc Avi kh L Az/oOr 1, 
LO 1 AyEL, AZ2ET, 
0 ; 

0. 0 i! Az 

| Rigs ee 1 


1 Vgl. FuBnote 3 auf S. 280. i 
2 Vel. S. 333 der in FuBnote | auf S. 280 zitierten Diss. 
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c) Die Verwindungsmatrix &%. Diese stellt eine Koordinatentransformation vom System 
(€,) auf das System (x, y) dar: 
COST) | sin y \ 


cos y 


8, =|— See = et 
—sin y cos y 
—sin y | cosy 
—sin y cos y 


sin | cos y 


Zur Riicktransformation dient die Kehrmatrix %71, bei der nur die Vorzeichen der Elemente 


sin y gegentiber ¥ vertauscht sind. 5 


d) Die Massenmatrix. Die Massenkrafte der Schwingung treten in den einzelnen Punkt- 
massen als jeweils neu hinzukommende Querlasten auf. Sie haben fiir die Masse i die GréBe 
uw S; Az y; fir die y-Richtung und yu w? S; Az x; fiir die x-Richtung mit yw als spezifischer Masse 
und S; als Querschnittsflache. Als Querlasten stehen die Werte yu jw? S; Az in der vierten und achten 
Zeile der Massenmatrix und in der ersten und fiinften Spalte, da sie den Verschiebungen y; und x; 
des Vektors (1) zuzuordnen sind. 

An den einzelnen Schaufelteilen greifen die Fliehkrafte C; an, welche die Schaufel in die Null-Lage 
zurickzubringen versuchen. Die Fliehkrafte stehen bei einer Auslenkung der Schaufel in y-Rich- 
tung parallel zur z-Achse (Abb. 3a), bei einer Auslenkung der Schaufel in x-Richtung dagegen in 
der radialen Richtung (Abb. 3b). Die Fliehkrafte werden fiir diesen Fall in zwei Komponenten 


zerlegt. Die Komponente C’ in z-Richtung ist der GréBe nach etwa gleich der Fliehkraft C, da die - | 


Schwingungsausschlage klein sein sollen. Die Komponente quer zur Biegelinie der Schaufel U 
liege, wie die Massenkrafte der Schwingung, in der Richtung der Koordinate x. 


Abb. 3. Fliehkraft an einem Schaufelteil — Abb. 4. Schaufel mit Fliehkriften in n Abschnitte unterteilt, — 
a) in der Ebene y; z; b) in der Ebene x; z. a) ganze Schaufel; b) Schaufel am Punkt i. 


In Abb. 4 ist die Schaufel in n Abschnitte mit dem Index » von 0 bis n unterteilt dargestellt. 
Die an den Massenpunkten angreifenden Fliehkrafte C, miissen von der obersten Masse m,, ange- 
fangen durch die darunter liegenden Abschnitte in den Schaufelfu8 geleitet werden. Im Abschnitt 
t+ 1/t sind alle Flichkrafte von C, bis C;,, aufzunehmen, im Abschnitt i/i—l alle Fliehkrafte 
von C,, bis C;. Mit der Rotationsgeschwindigkeit Q ist 


CG, & C= wus (RA) is (6) 
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Die innere Querkraft Q; ist im Stab proportional zu dex. zu ubertragenden Summe der Flieh- 
_ krafte und zum Neigungswinkel der Biegelinie gegen die z-Richtung. Die Komponente U, soll 
vorlaufig nicht beriicksichtigt werden. Die innere Querkraft ist unterhalb des Punktes i (Abb. 4b) 


QF pO ¥ S,(R +2) 4:(%) | 


v=n L 


und oberhalb des Punktes i (7) 


E oe dx) 
OF =p? » S,(R + z,) Az rae 


von 


Bildet man daraus den Mittelwert fiir den Ubergang tiber die Punktmasse, so erhalt man als An- 
derung des Momentes fiir die Lange Az 


1 ae | 
AM, = Qq:42 = |p OF S,(R+2)+ ¥'S,(R+2)|421(4).. (8) 
Der in der geschweiften Klammer stehende Ausdruck wird in der Massenmatrix in die siebte Zeile 
und sechste Spalte eingesetzt. Die entsprechende Gleichung gilt auch fiir die Auslenkung der 
Schaufel in y-Richtung : 


1 vues | dy 
AM, ;= {1 oe & S;(R+2) + > S,(R+ 2,)| z| a 


v=n 


; 0) 
Der Ausdruck in der geschweiften Klammer wird in die dritte Zeile und zweite Spalte der Massen- 

matrix eingesetzt. 
Nun ist noch die waagrechte Komponente der Fliehkrafte in x-Richtung zu bericksichtigen. 
Nach Abb. 4a besteht die Proportion 
: x; U. 


Ree Ce : (10) 
Dies ergibt fiir die Stelle i eine Querbelastung 
x. 1 

Ur= ae C; = {p 2? S, Az\ x. (11) 


Diese Querbelastung hat die gleiche Form wie die Massenkraft der Schwingung (es steht nur 
Q statt w) und die gleiche Richtung. Der Ausdruck in der geschweiften Klammer ist mit der Ver- 
schiebung x des Vektors (1) multipliziert. Er wird dem Element der achten Zeile und finften 
Spalte der Massenmatrix hinzugefiigt. 

Die Massenmatrix der rotierenden Schaufel lautet nun 


( 1 0 008 
0 1 C0: 
il 
pp $22 |S; (R + 3) 0 
O Boe 1 0 
4 ES, (R+2)] de 
“Sw? Az 0 Ort 
= eae ‘i il 0 0 0 
0 1 0 O 
: Gow 
; pO 5S, (R +2) 
2 0 pee 1 0 
| | ee er 
\ wS (w?+Q2) Az 0 Oead 
: (12) . 


e) Ubergang von der Stelle z = 0 auf die Stelle zs = L. Hierbei gilt, wenn man die ein- 
zelnen Abschnittsmatrizen (3) nach (2) aneinanderfiigt, _ 
be —MBLBYyn ar MBLBY 4M BLVHpi-M BEB Yo. (13) 
; - ; 19* 
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Man kann also durch Multiplizieren aller Abschnittsmatrizen die Matrizengleichung fiir die Be- 
ziehungen zwischen SchaufelfuB (z = 0) und dem Ende der Schaufel (2 = L) aufstellen mit der 
Gleichung (14), in welche die Randbedingungen eines einseitig eingespannten Stabes schon ein- 
gefiihrt sind. Dabei stellt die Matrix St das Produkt aller Abschnittsmatrizen dar: 


iat = Seo & (14) 


oder ausfihrlich geschrieben 


Y Ki Ky a re | Ss Ke a \ ( 0 \ 

dy /dz Koy Ky, Kez Kos | Kes Keg Kaz Keg 0 
0 K3, Kg, aKegg teaaa' Ky, Kg | Ky, Kgs! — M, 
: Ka Ken | Kas Kae) | Kas Kae | Bax Kas || |—@ re 
rye Ks, Keg ~ Kya Ky | Ks; Ks, Ks, Kes 0 

dx/dz Kot kss Kes Koyjpnildsee- tiga ie Kes 0 
0 Ky, Ky, | Keak Ky5 Ky, | K,, Ky, | —'M, 

0 Ky, Keg. Keg, Keg | Kz, Kg, | Ky, Kee |} \—Q. 


Daraus ist ein Gleichungssystem mit vier Gleichungen zu entnehmen, deren Koeffizienten ein- 
gerahmt sind. Die Bedingung fiir die Berechnung des Eigenwertes ist, wie bekannt, daB die Determi- 
nante der Koeffizienten verschwindet. Man rechnet also fiir mehrere Werte von w den Wert der 
Determinanten aus und sucht die Nullstellen auf. Die entsprechenden Werte von qm sind dann 
die Kigenfrequenzen. 


3. Die elastisch eingespannte Schaufel. Um die Elastizitat der Einspannung zu beriicksichtigen 
kann man am SchaufelfuB eine weitere Matrix einfiihren. Die Weichheit der Einspannung kann 
in Umfangs- und Achsrichtung verschieden sein. Es ist also im Punkt z = 0 eine Anderung der 
Neigung der Biegelinie, nicht aber eine Verschiecbung méglich. Die EinfluBzahl h sei die Anderung 
der Neigung unter der Wirkung eines Momentes. Damit ergeben sich fiir die einzelnen Kom- 
ponenten des Vektors (1) vor und hinter der Kinspannung 


1 =%=?9, x, =2 0; 
dy; __ dy as dx, dx 
rd aaa hy My,o» ieee M,.,0; (16) 
My,1 = Myo, Mi = Mi as 
Q,,1 = Qy,0 Op = Qu, 0» 
Die Matrix fiir den Schaufelfu8 lautet nun 
1 0 | 
if h 
ee 
1 0 
Q ie Weil 
= : 1 
od 0 (17) 
acne ae 
1 
Oi eel 


Diese Matrix steht in dem Matrizenprodukt nach Gl. (13) als erste Matrix am SchaufelfuB, also 
ganz rechts. 


4. Die Gelenkschaufel. In der Praxis wird eine in Umfangsrichtung gelenkige Befestigung der 
Schaufeln in gewissen Fallen angewendet. Dabei befinden sich in axialer Richtung angeordnete 
Gelenkbolzen in der Scheibe. 


a) Ideales Gelenk. Kin solches 1a8t sich durch Andern der Randbedingungen in die bisherige 
Rechnung mit und ohne Fliehkrafteinflu8 einfiihren.. In der x-Richtung sind am SchaufelfuB die 
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Verschiebung und das Moment null. Damit tritt an die Stelle der Matrizengleichung (15) die Glei- 
chung (18) als Beziehung zwischen dem Innen- und Au®enrand des Problems 


] Ky, Ky, Ky, Ky, |K; Ky, kK, sec ( 0 \ 
dy/dz Ky, Ke. mc Ky, | K,, Kz, Ky, Kye 0 

0 K;, Kg, | sew Kar ase thse) Key Ky | — M, 

0 Ky, Ky, | Kas Ky, | Ky; Kel Ky, |Kas| — Uy 

ue Ks, Ks. ss Kea | Ky, Ks, Ky Ks. ee - 
dx/dz Won Mey Ket pai. 3% Kp: Kyo Keg | | dxjdz 

0 Kj, Ky, Ky Kyq| | K,; Kye Ky, Kre | 0 

0 \Ke, Ky. | Kgs Key | Ky; Keo Kg, Kes] ms ) 

rete B 2==0 


Die eingerahmten Stellen bezeichnen die Koeffizienten des Gleichungssystems, aus dessen Deter- 
minante die Eigenwerte berechnet werden. 


b) Wirkliches Gelenk. Tatsachlich kénnen die Bedingungen eines Gelenks im héheren Dreh- 
zahlbereich nur durch eine Anordnung aufrechterhalten werden, welche reines Abrollen ohne Rei- 
bung gewahrleistet!. Bei einem Bolzen oder einem LavalfuB ohne Spiel wiirde die Reibung eine 
starre Einspannung erzeugen. Beim Abrollen verschiebt sich der 
Punkt der Auflage des Gelenks in Richtung x (Abb.5). Der 
Winkel ¢, welchen der durch den Berithrungspunkt gehende Radius 
des Bolzens 2 mit der Senkrechten bildet, kann durch den Bolzen- 
radius rz, durch den Radius des zu einer Schlaufe erweiterten 
SchaufelfuBes r, und durch die Neigung der Schaufel in «-Richtung 
folgendermaBen ausgedriickt werden: 


r, (dx 
e=p nla), 9° (19) 


Die Verschiebung e des Beriithrungspunktes ist damit 


ry, (dx 30 
ee la ee yal Ieee (20) Abb. 5. Gelenk. 


Das Biegemoment, welches durch das auBermittige Angreifen der Fliehkrafte der ganzen Schaufel 
mit Schaufelfu8 im Berithrungspunkt entsteht, ist 


M,,, =e ee (21) 
Setzt man die Ausdriicke fiir e und C in die obige Gleichung ein, so erhalt man 
at, ey dx 
Mya = {2% nor S 1S,(R+5) del (2) (2) 


Es 148t sich nun wieder eine Matrix des Gelenks aufstellen, indem man den Ausdruck in der 
geschweiften Klammer der Gleichung (22) der Neigung in x-Richtung zuordnet : 


1 
: 0 
1 
1 
x | (23) 

aie 1 0 0 0 

0 1 OpAg= 0 

0 rr. 
Op is S(RAE a Az ls 0 
0 0 Orne el 


1 Vgl. FuBnote 7 auf S. 280. 
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Diese Matrix wird in der gleichen Weise wie die Einspannungsmatrix (17) aus Ziff. 3 am SchaufelfuB, 
also in (13) ganz rechts, eingefiihrt. Nach dem Aufmultiplizieren gilt dann (18) fiir die Einfiihrung 
der Randbedingungen und die Aufstellung der Determinanten. Durch Variation von @ findet 
man wie in Ziff. 2 die Stellen von w, bei denen der Wert der Determinanten null ist. 


5. Die Schwingungsform. Die Berechnung der Schwingungsform zu den Eigenfrequenzen geht 
bei der eingespannten Schaufel von Gleichung (15) aus. Mit den fir eine Eigenschwingungszahl ra) 
errechneten Koeffizienten kann man fir die eingespannte Schaufel folgendes Gleichungssystem 


anschreiben: 


K35 a, + K5, he + Kz, M, + Kg, Q, = 0, 
Ky; M, hikga Q, + Ky, M, + Ky, Q, = 9, 
Ky M, a> K,4Q, nits Ky, M, a K,, Q. =0, 
Kgs M, ay Kgq Ue sik Kgy M, rile Kgg Q, = 0. 
Die Koeffizienten sind die gleichen, welche auch die Determinante zur Bestimmung von w enthalten 
hatte. Das obige Gleichungssystem ist homogen, so daB man die vier Unbekannten nur in ihren 


Verhaltnissen zueinander erhalt. Setzt man eine davon gleich 1, so benétigt man nur noch drei 
Gleichungen. Man kann die iibrigen drei Unbekannten aus folgendem Gleichungssystem bestimmen: 


(24) 


K33 M, + Ks, Q, + Ke, M, = — Ks, 
Ki, M, + Kis QO. + Bee M, = — Ky, Q,=1. (25) 
Ky M, + K,4Q, + Ky, M, =— K;,, 


Damit ist der Vektor fiir z = 0 bekannt. Man kann durch Multiplizieren mit der Matrix nach (15) 
den Vektor fiir z = L ausrechnen. Die Zwischenvektoren folgen aus den Produkten der ent- 
sprechenden Abschnittsmatrizen nach (13). Um die Schwingungsform zu erhalten, entnimmt man 
den Vektoren die Werte fiir y und x. Man kann in einem Koordinatensystem (y; x) die Schwin- 
gungsform als Kurve auftragen, wobei man die Stelle z an die Kurvenpunkte hinzuschreiben muB. 
Die Werte fiir y und x sind nur in ihrer relativen GréBe gegeben. 

Bei der Gelenkschaufel ist der Vorgang entsprechend (18) mit den dort eingerahmten Koeffi- 
zienten. Anstelle von — M, tritt hier dx/dz als Unbekannte im Vektor fiir z = 0 auf. Das dem 
Gleichungssystem (25) entsprechende System lautet fiir die Gelenkschaufel 


- 


= Sid a 

K33 M, te Ky4 Q, Tee K3 = ess K3. ? 

2 2 dx 

Kys My + Ky, Q, — Kye = — Kg, Q,=1. (26) 
dx 

Kj, My + Ky, Q,— Kye = — Kj,, 

Die vierte Gleichung 

dx = 

Kgs M, ci Kg, Oy Kge oo Kgs (27) 


wird auch hier nicht benétigt, kann aber verwendet werden, wenn man eine andere der obigen drei 
Gleichungen dafiir weglaBt. 


6. Zahlenbeispiel. Die Kigenfrequenzen einer Schaufel, welche in ihrer Biegsamkeit der Schaufel 
eines Axialverdichters der vorderen Stufen entspricht, wurden unter Annahme starrer Einspannung 
am Schaufelfu8 und unter der Annahme einer gelenkigen Befestigung nach dem beschriebenen ~ 
Verfahren berechnet. Die Daten sind in Tabelle 1 enthalten. 


% Tabelle 1. Daten der Schaufel 


t 0 Oia Sao 0,3 0,4 OSs jaac0y6 OF 4058 0,9 1,0 
S mm? 210 192 175 162 152 | 144 137 132 128 | 
Je mm‘ 713 517 392 307 244 197 165 140 122 
Jy mm4 30800) 28150 | 25700 | 23650 | 22000 | 20800 | 19820 | 19060 | 18500 | 18050 17700 
Ve Nee 88,7 85,5 80,0 74,5 69,5 64,8 60,0 99,6 51,8 48,0 45,0 


L= 250mm, 2Q— 628,27 sec 4, E = 0,72 - 104 kg/mm? 
Az= 25mm, R = 250mm, b= 2,854 - 19-20 “e/see 
mm 


Schw: 


oh th 
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GemaB Abb. 2 wird jeweils beim Abschnitt i—1/i die Flache von i eingesetzt und beim letzten 
Abschnitt die halbe Flache. Bei den Tragheitsmomenten und dem V erwindungswinkel wird das 
arithmetische Mittel zwischen i — 1 und i gebildet. Die so erhaltenen Werte stehen in Tabelle 2 


Ge 


Tabelle 2. Daten der Rechnung 


Abschnitt 0,.1>0,2 | 0,2+0,3 | 0,3 + 0,4 


S mm 1o2a3) nls) meena | 137 132 | 128 | 125 | 61,5 

Js mm! 615 | 454,5 | 349,5 | 275,5 | 220,5 | 181,0 | 152,5 | 131,0 116 106, 

Jn mm! 29475 | 26925 | 24675 | 22825 | 21400 | 20310 | 19440 | 18780 | 18275 | 14875 
° PO beets | 17,25. | T4saim orks, | 62,4 | - 57,8 53,7 | 49,9 46,5 


Ferner wurde die Rechnung mit der Annahme eines Gelenks im Schaufelfu® durchgefiihrt, um 
das Schwingungsverhalten einer Gelenkschaufel zu untersuchen, ohne auf die Méglichkeit und 
ZweckmaBigkeit in der Praxis einzugehen. Dabei waren die Radien r, = 13 mm und r, = 10 mm 
bzw. 12,5 mm. Um die gleichen Daten (Tabelle 1 und 2) verwenden zu kénnen, liege der Punkt P 
(Abb. 5) bei z = 0 der bisherigen Schaufel. 

Die Higenfrequenzen sind fiir die Schaufel mit starrer Einspannung und fiir die gelenkig be- 
festigte Schaufel in Tabelle 3 zusammengestellt. See 


Tabelle 3. Eigenfrequenzen in Hz 


starre Einspannung 


Rechnung 


Gelenk r, = 10 mm 


Gelenk r, = 12,5 mm 


Versuch 


starre 
Einspannung 


= 1/s 0 50 
Bee! | 78 * 105 
es: TF) 292. 325 
=~ III| 746 766 

IV| 959 993 


100 0 


156,5 0 
404 189 
816 670 

1089 852 


5 50 


TiO 
189,5 243.5 
67> 718 
852 876,5 


100 


134 
350,5 
807 
970 


0 
189 
670 
852 


5 


tS 
190 
671 
853 


50 100 


92,5 151 
281,5 388 
744 814 
915 1048 


0 


78 
286 


1020 


Fur die nicht rotierende, starr eingespannte Schaufel wurden die Eigenfrequenzen gemessen! 
und in der letzten Spalte der Tabelle 3 zum Vergleich zu den gerechneten Werten in der ersten 
Spalte angefiihrt. 


y-Richtung 


v 


20.012 
—0,051 
—0,119 
—0,218 
—0,335 
—0,497 
—0,668 
—0,852 
—1,043 
—1,236 


wv 


ooceo mos sssssss 


wie ss 
Okwhd 


x-Richtung 


rossss 
SOMOAD 


ey 


Tabelle 4. Verschiebungen y und x der nicht rotierenden Schaufel 


Starre Einsp: 


—0,0058 
—0,0211 
—0,0415 
—0,0605 
—0,0714 
—0,0687 
—0,0500 
—0,0172 

0,0250 

0,0709 


0,086 
0,348 
0,857 
1,705 
2,855 
3.990 
4,529 
3,951 


0,201 
0,560 
0,534 
—0,195 
—1,435 
—2,479 
—2,401 
—0,629 


2,658 
6,697 


1 P, Hermann, Konstruktion (im Druck). 


0,00027 
0,00106 
0,00208 
0,00252 
0,00161 
—0,00042 
—0,00171 
0,00027 
0,00650 
0,01539 


—0,00193 
—0,00614 
—0,01005 
—0,01083 
—0,00710 
—0,00032 

0,00571 

0,00728 

0,00355 
—-0,00318 


—0,001 
—0,003 
—0,010 
—0,027 
—0,059 
—H,102 
—0,190 
—0,290 


—0,082 
—0,154 
—0,207 
—0,230 
01 L 
—0,149 
—0,045 

0,092 

0,248 

0,411 


—0,0008 
—0,0037 
—0,0108 
—0,0240 
—0,0419 
—0,0573 
—0,0591 
—0,0381- 
0,0054 

0,0611 


—0,0239 
—0,0385 
—0,0355 
—0,0130 
0,0211 
0,0502 
0,0558 
0,0280 
—0,0287 
—0,0994 


annung Gelenkige Befestigung 
TED Iv II IIr Iv 


0,0004 
0,0015 
0,0030 
0,004.4 
0,0052 
0,0055 
0,0064 
0,0090 


—0,0050 
—0,0095 
—0,0120 
—0,0116 
—0,0079 
—0,0022 
0,0034 
0,0070 
0,0083 
0,0082 
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In Abb. 6 sind die Eigenfrequenzen iiber der Rotationsgeschwindigkeit Q/21% gee sale 
der eingespannten Schaufel sind die untersten vier Eigenfrequenzen als Kurven 1 4 ‘ au - 
zeichnet. Da es sich um biegsame Schaufeln handelt, ist der Einflu8 der Flichkraft bedeutend. 
Die Kurven gehen von der Ordinatenachse mit waagrechter Tangente aus. 


Hz AFigenfrequenz 
1100 g 7 


1089 
4 1048 
: l 
7000 F oa a a 
959 # eu 4970 
ZA 
if at ae 
ae 
900 V Gun “915 roe a 
a ap 57: 5 
ae Shae ae, 4 
652 “ 816 | ory 
800 ———a 807 
766 eae ee 
746 Es a 
I See he 
700 Se WS 
pee Te 
670 
600 
500 t 
404 
a 38 
FAROE 
305 se Been : 
Ep : = 
500 292. 2648 esa, Soe 
o a 
= i or 
aes) 
co Sa 
200 BES ye 
Ge 
189 1665 
157 
Nag, 
Ait Me | en 
100 E, _ = a 
78 le, 5 — 
- —l7 
é ES 
v6 = er 
4,91 Rotationsgeschwinaigkeit Blea 
5 a 
0 $5 50 100 1/8 
Abb, 6. Die Eigenfrequenzen der starr eingespannten und der gelenkig Abb, 7. Schwingungsformen der starr eingespannten Schaufel, — 
befestigten Schaufel. a) Grundschwingung, erste und zweite Oberschwingung; b) dritte 
Oberschwingung. 


Fir die gelenkig befestigte Schaufel sind die Kurven der Eigenfrequenzen fiir r, = 10 mm mit 
der Bezeichnung G; bis G, und fiir r, = 12,5 mm mit der Bezeichnung Gj bis Gj aufgetragen. Die 
Kurven der Gelenkschaufel nahern sich fiir wachsende Rotationsgeschwindigkeit den Kurven fir 
die eingespannte Schaufel. Am interessantesten ist dabei der Verlauf der Kurven G, und G/. Diese 
kommen nicht mehr mit allen erregenden Ordnungsstrahlen zum Schnitt. Die Kurve Gi, geht mit 
der Neigung 1,57, die Kurve Gj mit der Neigung 2,85 vom Nullpunkt aus. Die Kurve G, kommt 
also selbst mit der zweiten Ordnung der Drehzahl nicht mehr in Resonanz, die Kurve G{ nicht mehr 
mit der dritten Ordnung. Das giinstigere Schwingungsverhalten der Kurve G, mu aber mit einer 
gr6Beren Flachenpressung der beiden zylindrischen Kérper des Gelenks erkauft werden.1:2 


1 K. Léffler, VDI-Berichte 24 (1957) S. 27. 
2 G. Hirsch, Jahrbuch d. Wiss. Ges. f. Luftfahrt 1958, S. 174. 


\ 
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Nach Ziff. 5 sind die Schwingungsformen fiir Q = 0 mit Q, 9 = 1 fiir die eingespannte und die 
gelenkig hefestigte Schaufel berechnet worden. Die Verschiebungen y und x sind in der Tabelle 4 
enthalten. Bei der Gelenkschaufel ist die Kigenfrequenz der Grundschwingung null fiir Q = 0; 
damit entfallt es, die Schwingungsform I darzustellen. ae ; 


¥) 


=i 


Abb. 9. Schwingungsformen der starr eingespannten 
mit 2/2 «= 100 1/s rotierenden Schaufel. 


=4 ‘ 1 


Abb. 8. Schwingungsformen der gelenkig befestigten Schaufel. — 
a) erste; b) zweite; ¢) dritte Oberschwingung. 


Die Probe ist mit Hilfe der Bedingung der Orthogonalitat! 
w,.Gtwg=0 fir o+f6 


durchgefiihrt worden. Jede Spalte in der Tabelle 4 wurde dabei als Kigenvektor 1v mit zwanzig 
Komponenten aufgefaBt. Die Schwingungsgrade sind mit « und B bezeichnet. Die Massenmatrix © 


1 B. Jager, Ing.-Arch. 27 (1959) S. 51. 


Rey eye 
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ist durch den Verlauf der Flache S ither die Lange ¢ gegeben 


ty 

7,0 
ip 
Le 
7 
= if = 1 
Abb. 10. Schwingungsformen der gelenkig befestigten rotierenden Abb. 11. Schwingungsformen der gelenkig befestigten rotierenden 
° Schaufel (Q/2 x = 100 1/s; r, = 10 mm). Schaufel (2/2 2 = 100 1/s; r, = 12,5 mm). 


In Abb. 7 und 8 ist zu den Schwingungsformen die Koordinate z als € = 2/L hinzugeschrieben. 
Die Elemente der Kigenvektoren aus Tab. 4 sind zur Normierung jeweils durch die gréBte Kompo- 
nente dividiert worden. In Abb. 9, 10 und 11 sind die Schwingungsformen der rotierenden Schaufel 
aufgezeichnet. 


(Eingegangen am 6. November 1959.) 


Anschrift des Verfassers: Dr.-Ing. Berthold Jager, Stuttgart-Wangen, Kemptenerstr, 22 , 
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Analyse von Messungen turbulenter Grenzschichten 
mittels der Wandgesetze 


Von W. Szablewski 


1. Einleitung. Das logarithmische Geschwindigkeitsgesetz turbulenter Grenzschichten kann 
man bekanntlich aus der Bedingung t & t, (t Schubspannung, T, Wandschubspannung) auf der 
Grundlage der Prandtlschen Mischungsweghypothese erschlieBen. Die Bedingung besagt, daB 
innerhalb einer wandnahen Zone der Grenzschicht Tragheits- und Druckkrafte gegentiber den 
Reibungskraften vernachlassigbar sind. Die Vernachlassigbarkeit der Tragheitsglieder in der 
Bewegungsgleichung der Grenzschichtstrémung erscheint fiir eine wandnahe Zone aus der an der 
Wand zu erfiillenden Haftbedingung gerechtfertigt. Man erhalt dann zunachst durch Integration 
die Gleichung 


(1) 


(x,y Koordinaten langs und normal zur Wand; p statischer Druck). Die Existenzschicht des 
logarithmischen Gesetzes ist nun weiter gekennzeichnet durch Vernachlassigbarkeit des zweiten 


Gliedes in der Klammer gegen Eins. Man ersieht aber auch aus (1) unmittelbar, da8 fiir wachsenden 
dp/dx 


ran(t tity 


Parameter die Existenzschicht des logarithmischen Gesetzes immer weiter eingeengt und 


0 
dieses schlieBlich ausgeléscht wird. Damit ist evident, da8 im Bereich turbulenter-Grenzschicht- 
strémungen mit starkem Druckanstieg und kleiner Wandschubspannung ein anderes als das log. 
Wandgesetz zu erwarten ist. 

Mit Ausgang von Gleichung (1) und auf der Grundlage der Prandtlschen Mischungsweghypo- 
these habe ich in zwei Arbeiten! ein dem EinfluB des Druckgradienten Rechnung tragendes ,, Wand- 
gesetz turbulenter Grenzschichten mit Druckanstieg“ aufgestellt. Der Vergleich mit den von 
J. Nikuradse? durchgefihrten Messungen turbulenter Strémungen in divergenten Kandlen ergab 
gute Ubereinstimmung. In der vorliegenden Arbeit wird der Vergleich mit neueren Messungen 
von Profilstrémungen durchgefiihrt, die gleichfalls eine gute Bestatigung des genannten Wand- 
gesetzes ergeben. 

Einige prinzipielle Bemerkungen zur Prandtlschen Hypothese erscheinen an dieser Stelle 
angebracht. Ihren ersten Erfolg zeitigte diese Hypothese in der Herleitung des log. Gesetzes.? Die 
Existenz dieses Gesetzes kann man allerdings auch hypothesenfrei erschlieBen, wenn man an- 
nimmt, daB einerseits ein von den auBeren Bedingungen unabhangiges Wandgesetz und anderer- 
seits ein von der Zahigkeit unabhangiges Mittengesetz existiert*. 

Doch spricht dieser Umstand natiirlich nicht gegen die Prandtlsche Annahme. Da® man aus 
ihr das logarithmische Gesetz ableiten kann, stellt eine Stiitze der Hypothese dar. Wenn nun 
_ weiter die vorliegende Untersuchung eine Bestatigung des aus (1) auf der Prandtlschen Grundlage 
_ abgeleiteten Wandgesetzes fiir turbulente Grenzschichten mit. Druckanstieg ergibt, so muB man 
darin eine Konsolidierung der Prandtlschen Hypothese, zumindest fiir die wandnahe Zone turbu- 
lenter Grenzschichten, erblicken. : 

Die Prandtische Formulierung ir die turbulente Schubspannung 


au\2 
Tturb. = 0 B (5) (2) 


(o Dichte, u Geschwindigkeitskomponente in x-Richtung) mit dem Mischungsweg ! beruht letztlich 
auf der Lokalitatshypothese, daB namlich die turbulente Schubspannung bestimmt wird durch 
den lokalen Verlauf der mittleren Geschwindigkeitsverteilung u(y), d.h. durch® du/dy. Gegen 
diese Annahme der lokalen Bedingtheit sind mehrfach Bedenken, insbes. in der angelsachsischen 


1 W. Szablewski, Ing.-Arch. 22 (1954) S. 268, zitiert als Arbeit I; Ing.-Arch. 23 (1955) S, 295, zitiert als 
Avbeit II. | : 
" J. Nikuradse, Untersuchungen iiber die Strémung des Wassers in konvergenten und divergenten Kandlen, _ 
VD1.Forschungsheft 289 (1929). 
3}. Prandtl, Z. angew. Math. Mech. 5 (1925) S. 137. ; 
4 C, B. Millikan, V. Intern. Kongr. f. Mechanik, Cambridge, Mass. 1938. 
5 Vel. z. B. W. Szablewski, Ing.-Arch, 20 (1952) S. 67, 
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Literatur, erhoben worden. Was nun Phdnomene der sog. freien Turbulenz (Diisenstrahlen und 
Nachlaufstrémungen) mit ihren verhaltnismaBig groBen Mischungswegen betrifft, so ist dort 
allerdings, wie L. Prandtl' selbst bemerkte, die der lokalen Bedingtheit entsprechende differentielle 
Formulierung des turbulenten Austausches A, entspr. dem Ansatz tT = A du/dy, nicht als sachgemab 
anzusehen; Prandtl hat dort konstanten Austausch iiber dem Querschnitt vorgeschlagen. Ahnlich 
kénnten nach den Untersuchungen von F. H. Clauser® die Dinge hinsichtlich der auSeren Zone 
turbulenter Grenzschichten liegen. Fiir die wandnahe Zone turbulenter Grenzschichten jedoch ist 
die Hypothese der lokalen Bedingtheit der turbulenten Schubspannung und damit der Prandtlsche 
Ansatz (2) nach den bisherigen Ergebnissen als eine sachgemaBe Grundlage anzusehen. 

Eine andere prinzipielle Bemerkung betrifft die empirischen Koeffizienten der Wandgesetze 
[vgl. Abschn. 2, Gleichung (3)], die mit denen des log. Gesetzes identisch sind. Man findet haufig® die 
Ansicht vertreten, daB beihydraulisch glatter Wand diese K oeffizienten x und C bzw.v,,69/y universelle 
Konstanten waren. Da®B man in der Literatur schwankende numerische Angaben fiir diese Koeffi- 
zienten findet, braucht dem zunachst nicht zu widersprechen in Anbetracht dessen, dafs Mef- 
schwankungen nicht auszuschlieBen sind und daB die Auswertung fiir die Koeffizientenbestimmung 
keine eindeutige Prozedur ist. Es gibt jedoch zu bedenken, daB es auch Werte gibt, die ganz aus 
dem Rahmen der iiblichen Werte fallen; z. B. ergab die Auswertung* der Nikuradseschen® Messung 
der Strémung des Wassers in konvergenten Kanilen, die als besonders stérungsfrei anzusehen ist, 
die Werte x = 0,6; C = 9,3 baw. v,69/y = 9,5. Wir glauben daher, daB x und C bzw. v,0,/y eher 
als Modellkonstanten anzusehen sind, die von der Intensitat der Freistromturbulenz 
abhangig sind. Diesbeziigliche systematische Untersuchungen fehlen allerdings bisher. Handelt 
es sich um Strémungen an rauhen Wanden, so ist des weiteren die starke Abhangigkeit des Koeffi- 
zienten C baw. v,,69/y von der Wandrauhigkeit zu beachten®. 

Zum Vergleich mit der Theorie (Abschn. 3) haben wir die Messungen von A. Fage’ und V. A. 
Sandborn and R. J. Slogar® herangezogen. Unter der Voraussetzung, da diese Messungen ein- 
wandfrei sind, sind wir dabei so vorgegangen, daB wir unter Benutzung der experimentellen 
T Werte durch Anpassen der theoretischen an die experimentellen Geschwindigkeitsverteilungen 
die Koeffizienten x und C bzw. v,d,/y fiir das jeweilige Modell bestimmten. Diese erwiesen sich im 
Rahmen kleiner unregelmaBiger Schwankungen fiir das jeweilige Modell als Konstanten. Mit 
diesen Konstanten ergibt sich dann gute Ubereinstimmung der Theorie mit dem Experiment 
sowohl hinsichtlich der Geschwindigkeitsverteilungen als auch hinsichtlich der Wandschubspan- 
nungen Tp. 

In Abschn. 4 haben wir schlieBlich die bekannte Messung von G. B. Schubauer and P.S. Klebanoff® 
herangezogen. Die in diesem Report angegebenen Werte der Wandschubspannung t, sind um- 
stritten™. Fiir hydraulisch glatte Wand erscheinen sie als zu groB; es bleibt noch die Alternative, 
daB die Wand eine gewisse Rauhigkeit aufwies. Die vom Standpunkt der im vorliegenden Bericht 
gegebenen Theorie aus lésbar erscheinende Aufgabe, diese Frage zu entscheiden, erwies sich in der 
praktischen Durchfiihrung infolge der Empfindlichkeit der Rechnung und im Zusammenhang mit 
der Streuung der Messung als schwierig. Das erhaltene Resultat, daB die im Report angegebenen 
T- Werte zu groB sind, ist als nur bedingt richtig anzusehen. 


2. Die Wandgesetze. Das in Arbeit I abgeleitete Wandgesetz fiir turbulente Grenzschichten 
mit Druckanstieg lautet 


220g fee apie 4y —1 v, v, 6 | 
w= [2(Vi + = y—1)+m40—4], Be FInd + 22%) (3) 
mit der Schubspannungsgeschwindigkeit v, = /t/o und dem Formparameter 
_, » dp/dx | 
Se el ae 


1 L. Prandtl, Z. angew. Math. Mech. 22 (1942) S. 241, 
ee C. ee erence in Appl. Mech. IV (1956) S. 1. 
gl. z.B. D.A. Spence, J. Aeron. Sci. 23 (1956) S.3. Spence gibt an x = 0,4 = = 
4 W. Szablewski, Ing.-Arch. 20 (1952) S. Ay : Oe oie batiel eh aice ae ee ee 
5 J. Nikuradse, a. a. O. 
6 Vel. J. Rotta, Ing.-Arch. 18 (1950) S. 277. 
? A. Fage, R. &. M. 1852 (1938). 
° V. A, Sandborn and R. J. Slogar, NACA TN 3264 (1955). 
* G. B. Schubauer and P. S. Klebanoff, NACA Rep. 1030 (1951). 
10 Vel. z. B. J. Rotta, Mitt. aus dem Max-Planck-Inst. f. Strémungsforschung Nr. 8 (1953). 
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dy ist ein MaB fiir die Dicke der laminaren Unterschicht (0 Dichte, uw Zahigkeit, y = i/o kinematische 
Zahigkeit). 


In dimensionsloser Form kénnen wir auch schreiben 


Tete re ee REY, 1 4y —1 1 Pe,’ 
ale Rae 1)4 Pe ae a Lpln dn $e), (4) 
Fir F'— 0 ergibt sich aus (4) wieder das log. Gesetz 
ye 1 by 1 v, O 
in 4 al In ie tu). (6) 


Das Gesetz (3) wurde aus (1) unter Zugrundelegung der Prandtlschen Formel (2) abgeleitet, 
wobei in Modifikation des iiblichen Ansatzes der Existenz der laminaren Unterschicht durch den 
Ansatz | = x (y — 6,) Rechnung getragen wurde!. Man erhalt dann fiir das log. Gesetz die in (5) 
angegebene Form mit einer analysierten Fassung der Integrationskonstanten 


C=1(—1 pind tx), (6) 


% 


Doch ist es fiir die vorliegende Untersuchung belanglos, ob man C oder v, 69/y als empirischen 
Koeffizienten betrachtet. 

Gesetz (3) wurde unter der Voraussetzung abgeleitet, da® in der Strémung eine Schicht mit 
log. Gesetz zumindest noch rudimentaér vorhanden ist. In Ablésungsnahe (Ablésung v, = 0) 
ist das angegebene Gesetz (4), wie in Arbeit IT dargelegt wurde, im wesentlichen durch Einfithrung 
der Schubspannungsgeschwindigkeit 

Vi =| aes 
und demgemaf Einfiihrung der Reynoldsschen Zahl V, 6,/y der laminaren Unterschicht zu modifi- 
zieren. Wir gehen hier nicht weiter darauf ein, wie wir iberhaupt das Verhalten turbulenter Grenz- 
schichten in Ablésungsnahe ein anderes Mal untersuchen wollen. 


3. Vergleich mit Messungen?. a) A. Fage, Profile and skin-friction aerofoil drags, 
R. & M. 1852 (1938). Fiir die Anwendung der Formel (3) benédtigen wir auBer t) noch den Druck- 
gradienten dp/dx, den wir durch graphisches Differenzieren der experimentellen p-Verteilung 
erhielten (Abb. 1). 

Indem wir von der Vorstellung (vgl. Abschnitt 1) ausgehen, daB die empirischen Koeffizienten x 
und v,, d9/v (bzw. C nach Formel (6)) keine — aus anderen Messungen bekannten — universellen 
Konstanten, sondern Modellkonstanten darstellen, haben wir zundchst diese Koeffizienten be- 
stimmt. Das geschieht hier zweckmaBigerweise so, daB die experimentellen Geschwindigkeits- 
verteilungen u bzw. y = u/u, (u, Geschwindigkeit der Potentialstrémung) iiber der Koordinate 


Ss dp] ee 
214+ 2% y_1) +S (7) 
aufgetragen werden, wobei wir t) der Messung entnehmen; denn iiber & aufgetragen erscheint 
unsere Funktion (3) als Gerade! Man hat dann nur eine Gerade durch die wandnahen MeBpunkte 
au legen (vgl. Abb. 4), wobei allerdings die sehr wandnahen Punkte der laminaren Unterschicht 
oder der Ubergangsschicht, die als solche leicht er- anja 
kennbar sind, auszuschlieBen sind. Aus den Koeffizien 7 te 
ten der Geraden lassen sich dann nach (3) x und v, do/v 
ermitteln. Auf diese Weise erhielten wir fiir die einzelnen gj; 
Schnitte der Messung Re = 1,665- 10 von A. Fage 


A.Fage, R. & M. 1852 (1938) 


0,10\- 


l Castel tS 
10 4a 30 Go Jott 


i W, Szablewski, Z. angew. Math. Mech. 31 (1951) S. 131; ferner J. Rotta, a. a. O., 8. 3. a 
2 Rechnungen und Zolcbuitingeh der Ziff. 3 und 4 hat mein Mitarbeiter, Herr Herbert Moch, durchgefihrt. 
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4 Cr also eine Reihe von Werten mit echter Schwan- 
G005 | kung um die Mittelwerte 
v, Oo 

f x = 0,48 hess 6,71 (bzw. C =29,99) 5 
0,004 8 a 

8 3 Wir haben damit gleichzeitig den Nachweis er- 

. bracht, daB die Koeffizienten x, v, dg/v (bzw. C) 

ead : Konstanten des betrachteten Modells sind. (Die 

© 


obigen Schwankungen diirften im wesentlichen 
dadurch bedingt sein, da® die geschilderte Er- 
» Experiment A.Fage, R.& M. 1852 (1998) mittlung der Koeffizienten keine eindeutige Pro- 
Tie zedur ist). 
gor Nachdem einmal die beiden empirischen Ko- 
ie effizienten des Modells ermittelt worden sind, 
¢ ergibt sich gute Ubereinstimmung der Theorie 
OY 05 G O7 08 09 70° mit dem Experiment sowohl hinsichtlich der Ge- 
schwindigkeitsverteilungen als auch hinsichtlich 
der Wandschubspannungen Tp. 

Die Ermittlung der Wandschubspannungen t) aus den experimentellen Geschwindigkeits- 
verteilungen erfolgt dann in einfacher und sicherer Weise durch Anpassen der theoretischen Ver- 


G002 |- 


Abb. 2. 
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ers og. Geserz Pere 
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Theorie Cp=0,00353 
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Abb, 4. 


teilungen (3) an die iiber In y aufgetragenen experimentellen Verteilungen (vgl. Abb. 3). Die so 
theoretisch ermittelten Werte 1, stehen in guter Ubereinstimmung mit den experimentellen Werten; 
? 


vgl. Abb. 2, wo wir cy = To/ (5) uj aufgetragen haben. 
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Zu bemerken ist noch, da®B es demnach bei bekannter statischer Druckverteilung prinzipiell 
modglich ist, im Druckanstiegbereich aus der experimentellen Geschwindigkeitsverteilung die 
Wandschubspannungen 1) mittels Formel (3) zu gewinnen, sofern einmal die Modellkonstanten 
%, V, 0o/v (bzw. C) bekannt sind. Die Ermittlung derselben erfolgt in einem praktisch vorliegenden 


Fall am besten wohl an Hand von im Druckabfallbereich durchgefiihrten Messungen (einschlieB- 
lich T)) mittels des log. Gesetzes (5). 


Inwiefern nun das Wandgesetz (3) bzw. (4) die experimentellen Geschwindigkeitsverteilungen 
in einer wandnahen vollturbulenten Zone der Grenzschicht beschreiben, zeigen die Abb. 3 und 4. 
Dabei wurde in die Formeln mit den theoretischen t, -Werten elngegangen. 


b) V.A.Sandborn and R.J.Slogar, Study of the momentum distribution of 
turbulent boundary layer in adverse pressure gradients, NACA TN 3264 (1955). 
Der Vergleich von Theorie und Experiment erfolgt nach dem in a) geschilderten Vorgang. Abb. 5 
zeigt die durch numerische Differentiation erhaltene Verteilung des aueren Druckgradienten. 


4, ap/ax 
Ft pus 1 
004 '-, 
= ‘SIV = 
° 
Sandborn and Slogar, g003+ ° 
NACA IN 3264 (7955) 
G05 ie @) SU 
1S 002 + 
t © Experiment Sandborn and Slogar, Cae 
NAGA TN 3264 (1955) 


+ 0001 om 
° /Neorle 
Sian, 
2 0 i l \ { L 


5 6 7 é g ott 70 60 90 100 1) 120 790 inch 


Abb. 5. Abb. 6, 


Aus den iiber der Variablen € (Formel (7)) unter Benutzung der experimentellen t,)-Werte! 
aufgetragenen experimentellen Geschwindigkeitsverteilungen ergaben sich als Koeffizienten der 
durch die wandnahen Punkte gelegten Geraden fiir die MeBschnitte (vgl. Abb. 8) 


Vy Oo 


mit den Mittelwerten 
x =0,41 und ms == 5.0 ' (baw. CG =.3577) 
als Modellkonstanten. 


Nachdem einmal die Modellkonstanten bekannt sind, ergeben sich mittels Durchlegen der 
theoretischen Verteilungen (3) durch die experimentellen Verteilungen (vgl. Abb. 7) die in Abb. 6 
dargestellten theoretischen Werte der Wandschubspannung cy im Vergleich mit den experimentellen. 


Den Vergleich der Geschwindigkeitsverteilungen zeigen Abb. 7 und 8, wobei wieder in die 
Formeln mit den theoretischen Werten tT, eingegangen wurde. 


1 Entnommen der im TN angegebenen MeSkurve. 
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Anmerkung: Man kann dem Wandgesetz (3) auch die universelle Formulierung 


5) 
Se hart Vin 4 eee ‘| | 


Oh x x v 

mit | 

—- = (8) 
dp/|dx LD PRs | 

Sie) SS —1 + ln — ——— 
pe aa(/i+ 2A y— +n paz 
Ee * nach Formel (7 
é*=€+In “7 nach Forme (7) 

geben. 

ie 
40\- 
a Messung Sandborn and Slogan « = Messung Sandborn and Slogan 
NACA TN 326% (1955) NACA TN 3264 (71955) 
Theorie Cp=0,00245 Theorie Cp=G00746 
——— log. Gese/z ——— log. Gesetz r Ye 
mS . “ az rad ia 
SA oer 570 je eee See es 
i 4 m . Pee Eee 
rae Station II ae 
, L ae ; Station IV 
oe | pre 4 E 2 
| i J ul l [is \ | 
= =) Wa A 0 7 2 = oo. =2 =f 0 if é é 
Abb. 7. 


© Station I 


I Hea Sandborn and Slogan, 
f Theorie 4 It | NACA TN 326% (1955) 


——— log. Gesetz 


In der demgemaBen Auftragung fallen dann die Geraden der Abb. 4 und 8 in jeweilig eine 
Gerade zusammen und entsprechend werden die MeBpunkte gesammelt. Inwieweit das Wandgesetz 
fiir den einzelnen Schnitt erfillt ist, ersieht man jedoch besser aus den Auftragungen der Abb. 4 
und 8. Wir haben deshalb auf die Darstellung der universellen Formel hier verzichtet. 


4, Analyse der Messung von G. B. Schubauer and P. S. Klebanoff ,,Investigation of separation 
of the turbulent boundary layer“. NACA Rep. 1030 (1951). Wir stellten uns (vgl. Abschnitt 1) die 
Aufgabe, die Alternative zu entscheiden, ob die im Report angegebenen Werte der Wandschub- 
spannung T, als zu gro oder als durch Wandrauhigkeit bedingt anzusehen seien. 
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Es galt hier also, aus den wandnahen experimentellen Geschwindigkeitsverteilungen allein 
(ohne Zuhilfenahme wenigstens einiger gemessener t,-Werte zum Zwecke der Bestimmung der 
Modellkonstanten, vgl. Abschn. 3) die Wandschubspannungen zu ermitteln. 

Vom theoretischen Standpunkt aus erscheint die Aufgabe als sehr einfach: Bildet man im 
Giltigkeitsbereich des wandnahen Gesetzes (3) fir drei Punkte I, II, III eines Schnittes mit den 
Koordinaten yy, yr, yr z. B. das Verhaltnis 


ao ; 1355 8 
4 Some se St SES (€ nach Formel (7)), 
Yr — Pir /theor. fy — Sr 
so sind damit formelmaBig die empirischen Koeffizienten x, v, 69/v (baw. C) eliminiert und das 
Verhaltnis enthalt (bei bekanntem Druckgradienten dp/dx) nur die Unbekannte 7. 
Man hatte dann nur noch aus der experimentellen Geschwindigkeitsverteilung des Schnittes 
die Werte 9, Vy Py fiir die drei Koordinaten y,, Yip Nyy abzulesen; und aus der Gleichung 


ei Tinea UE ss Sy — Sur (9) 
Pr — Pr Jexp. = §¢ — Sat 


kénnte dann das gesuchte t, berechnet werden. Dabei wiirde im weiteren Verlauf der Rechnung, 
siehe Formel (12), sich Wandrauhigkeit an kleinen numerischen Werten v,, J9/v bzw. C zeigen (vgl. 


J. Rottat.). 
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205 t) C4 
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Abb. 9. 


In der praktischen Durchfihrung ist das leider nicht so einfach. Ks zeigt sich namlich, daB eine 
sichere Bestimmung der 1,-Werte aus (9) eigentlich ,,mathematisch genaue~ experimentelle Werte 
Pr Pir Pri erfordern wiirde; denn kleine Anderungen der y-Werte bedingen bereits groBe Ande- 


_rungen der 1)-Werte. Es stellte sich schlieBlich nur folgender Weg als gangbar heraus: Zunachst 


A s ; 
einmal galt es, durch die mit den normalen Schwankungen behafteten MeBpunkte eine ,,mittlere* 


1 J. Rotta, a. a. O., S. 3. 
20 


pel gs 


—, 
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MeBkurve zu legen. Zu dem Zweck erwies es sich als vorteilhaft, die MeBpunkte iiher der Koordi- 
nate /y aufzutragen. Diese Auftragung wird nahegelegt, wenn man beachtet, da nach Formel (3) 


sich ergibt 
w 2 lV (el) a peal 
one VCE) — sl Wow! 
34 Inde ting ta! fir PM (10)! 


0 


iv 
x Uy 


-- 


der Einflu8 der Reibung beginnt bei der Auftragung tber Vy, wie hier noch angemerkt werden 
soll, dort, wo die MeBkurve anfangt, sich von der Geraden abzukriimmen. 

Aus dieser, den wandnahen Bereich auseinanderziehenden, Auftragung (Abb. 9) lieBen sich 
nun die aus dem Trend der Messung fallenden Punkte in Wandnahe besonders gut erkennen und 
eine mittlere MeBkurve ermitteln. (Bei der Auftragung iiber In y sind die herausfallenden Punkte 
nicht so gut erkennbar). 

Die weitere Untersuchung erstreckte sich dann auf die vier Schnitte x = 19, 20, 21, 22 ft. Als 
Punkte J, IJ, III wurden fiir diese Schnitte die Punkte mit den Abszissen Vy; = 0,25, V¥_ = 0,11, 
Vn = 0,1 ft festgelegt, die als innerhalb des Geltungsbereiches des Wandgesetzes (3) gelegen anzu- 
sehen sind; vgl. Abb. 13. Bei diesen Abszissen wurden auf den jeweiligen MeBkurven der Abb. 9 
die Ordinaten@,, Y1;, Py; abgelesen und unter Mittelung die Kurven 9%, ;:Py, uber x aufgetragen; 
Zahlenwerte gibt die Tabelle. 


Pmittel 


x[fe] 

19 :  0,6100 
19,5 0,788 0,699 0,596 0,7884 | 0,6982 0,5916 
20 0,771 0,685 0,582 0,7672 | 0,6762 | 0,5708 
20,5 0,735 0,648 | 0,542 0,7414 | 0,6500  0,5460 
21 0,710 0,620 0,510 0,7108 | 0,6208 | 0,5176 
21,5 0,677 0,592 0,487 0,6766 | 0,5880 | 0,4874 
22 0,638 0,548 | 0,442 0,6388 | 0,5514 | 0,4536 
22,5 0,602 | 0,523 0,422 0,5974 | 0.5116 0,4166 
23 0,550 0,470 | 0,373 0,5500 | 0,4686 0,3754 
23,5 0,509 0,430 0,340 0,4954 | 0.4186 0,3296 
24 0,421 | 0,354 0,283 |. 0,4316 | 0,3600 | 0,2800 


Es wurden dann unter nochmaliger geringfiigiger Glattung die Kurven ¢,, — ~,,, sowie 9; — Pyy 
aufgetragen; bei Ablesung von den letzteren erhielt man 


lc et 


x| ft] iy Put 


Me Pr 
19 0.1077 | ° 0,1969 0.5470 
20 0.1058 _ 0,1962 0,5395 
21 0,1027 | 00,1926 0,5335 
22 0.0982 |  0,1856 - 0,5291 


Bei Einsetzen des ermittelten dp/dx (Abb. 11) ergab 
dann die graphische Liésung der Gl. (9) gemaB Abb. 10 


fiir -v,/u, = / too uj die Werte : 


Dy, E x[ jt] * 
ae 
G00 G0e5 G30 19 0,02850 
20 0,02570 
Abb. 10. 21 0,02377 
22 0,02256 


Aus der letzteren Abbildung ist die Empfindlichkeit der Rechnung in bezug auf kleine Anderungen 
der experimentellen y-Werte schon ersichtlich. Sie wird besonders deutlich, wenn man unter 


: Fiir die Ablosung wiirde aus (10) folgen: p = = V() Vy. Diese Formel hat jungst B.S. Stratford, J. 
Fluid Mech. 5 (1959) S. 1, seiner Betrachtung der Ablésung turbulenter Grenzschichten zu Grunde gelegt. 
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Festhaltung unserer y,- und Py- Werte sich die y,,;-Werte berechnet, mit denen (9) die im Report 
angegebenen wesentlich groSeren t)-Werte (vgl. Abb. 12) liefern wiirde. Es ergibt sich als 
»Alternativ-Differenz im Mittel 0,008 in Richtung kleinerer Py," Werte. Wenn ein um diese 
kleine Differenz von ~ 0,01 niedriger gelegener Verlauf von Py; auch nicht aus den Messungen 
ohne weiteres ablesbar erscheint, so ist sie andererseits doch — wie ein Blick auf Abb. 9 und die 
Tabellen zeigt — zu gering, um unser Resultat als sicher anzusehen. Es kann ihm daher nur ein 
bedingter Wert zugesprochen werden. 
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wobei das einzusetzende dp/dx durch numerische Differentiation gewonnen wurde (Abb. 11), 
ergab sich dann 


19 0.4105 8,374 
20 0,4055 9,230 
21 0,4152 9,160 
22 0,4408 8,102 
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mit den Mittelwerten 


~=042 und “8,72 (baw. C=7,51). 


Unter Zugrundelegung dieser Werte fiir x, v, 09/y als Modellkonstanten haben wir dann nach dem 
in Abschn. 3 dargelegten Verfahren den Vergleich von Theorie und Experiment auch fir dieses 
Modell durchgefiihrt. Abb. 12 zeigt die theoretischen und die im Report angegebenen Wandschub- 
spannungen. Abb. 13 und Abb. 14 enthalten die Geschwindigkeitsverteilungen. 
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Das Herausfallen der sehr wandnahen Punkte in den letzteren Abb. diirfte dabei wohl, wie 
man aus Abb. 9 ersehen kann, zu Lasten des Experiments gehen. 


5. Zusammenfassung. Es werden fiir turbulente Grenzschichten mit Druckanstieg experi- 
mentelle Geschwindigkeitsverteilungen in Wandnahe verglichen mit dem Wandgesetz dieser Ver- 
teilungen. Es ergibt sich gute Ubereinstimmung. 


(Eingegangen am 5. Dezember 1959) 


ae Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, Mohrensitr. 39. 
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